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Abstract One of the challenges of the modern photonics is 
to develop all-optical devices enabling increased speed and 
energy efficiency for transmitting and processing information 
on an optical chip. It is believed that the recently suggested 
Parity-Time (PT) symmetric photonic systems with alternating 
regions of gain and loss can bring novel functionalities. In such 
systems, losses are as important as gain and, depending on 
the structural parameters, gain compensates losses. Generally, 
PT systems demonstrate nontrivial non-conservative wave in¬ 
teractions and phase transitions, which can be employed for 
signal filtering and switching, opening new prospects for active 
control of light. In this review, we discuss a broad range of prob¬ 
lems involving nonlinear PT-symmetric photonic systems with 
an intensity-dependent refractive index. Nonlinearity in such PT 
symmetric systems provides a basis for many effects such as the 
formation of localized modes, nonlinearly-induced PT-symmetry 
breaking, and all-optical switching. Nonlinear PT-symmetric sys¬ 
tems can serve as powerful building blocks for the development 
of novel photonic devices targeting an active light control. 
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1. Introduction 

The concept of Parity-Time (PT) symmetry first emerged 
in the studies of quantum operators and their spectra [1-3]. 
Whereas it is well known that Hermitian operators have 
real spectra, Carl Bender and collaborators posed a ques¬ 
tion whether there exist non-Hermitian Hamiltonians which 
still have real spectra. They found that this can happen for 
a broad class of Hamiltonians which are symmetric with 
respect to the so-called PT transformations, defined by the 
parity operator P as P\j/(r,t) = i/a(— r,f), and time rever¬ 
sal operator T, as Ti//'(r,/ : ) = yC( r , — t), where \j/(r,t) is a 
wave function in quantum mechanics. When a simple one¬ 
dimensional system is described by the Schrodinger-type 
Hamiltonian with a complex potential U (x ), 

d 2 

# = --^ 2 +£/(*), U(x)=V(x) + M(x) (1) 

where V (x) and W (x) are both real, the system (1) is called 
PT-symmetric provided U*(x) = U (— jc), i.e. V(x) = V(—x) 
and W(x) = —W(—x). Although initially the analysis was 


motivated by the development of rather abstract quantum 
theories [1-3], it soon emerged that the concept of PT sym¬ 
metry is relevant to various physical systems. 

In optics, the propagation of light can often be described 
theoretically as an evolution governed by an effective Hamil¬ 
tonian (1) in the so-called paraxial approximation. In this 
case variation of the medium refractive index is described 
by the potential U (x ). Hermitian Hamiltonians correspond 
to the cases where the optical energy is conserved and U (x) 
is real. On the other hand, the presence of losses or gain in 
optical structures maps to non-conservative non-Hermitian 
operators (1) with a complex U(x). However, if such an 
operator possesses PT-symmetry and has a real spectrum, 
then the amplitudes of optical modes would be preserved, 
corresponding to an effective compensation of the effects 
of loss and gain, which can be important for a range of 
practical applications. First theoretical proposals for optical 
PT systems were formulated a decade ago [4,5], and it was 
shown that in waveguiding structures, P operator performs 
to spatial reflection, whereas T reverses the propagation 
direction. Accordingly, conventional optical PT-symmetric 
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structures have symmetrically positions regions of loss and 
gain, corresponding to the same real but opposite imaginary 
parts of the optical refractive index. The PT symmetry defini¬ 
tion can be generalized to include the gauge transformations 
enabling applications to practical structures without exact 
balance of gain and loss [6], and including the structures 
with pure loss [7]. 


The spectrum of PT operators is not always real, and is 
inherently connected to the symmetry of the eigenmodes. 
If an eigenmode profile is invariant with respect to the PT 
transformation, then it has real eigenvalue. However, the 
modes with PT-symmetry-broken profiles generally have 
complex eigenvalues. The PT-symmetry breaking transition 
and the associated change from real to complex eigenvalues 
can be generally observed when the amount of gain/loss is 
increased [6,7]. Such a phase transition can in particular 
facilitate single-mode lasing [8,9] and find various other 
applications (see the recent review paper [10] and references 
therein). 


Even before the PT symmetry concept was introduced, 
it was suggested that optical systems with gain and loss can 
improve the operation of nonlinear all-optical switching de¬ 
vices, including the lowering of the switching power [11,12]. 
The development of PT concept brought a new perspective 
and stimulated further active investigations of nonlinear ef¬ 
fects in PT-symmetric systems. In addition to switching, 
the topic of nonlinear self-focusing and the formation of 
self-localized solitons received a strong attention. Optical 
solitons in the conservative systems or the systems with low 
losses originate from a balance of nonlinearity and disper¬ 
sion and usually form families of localized solutions [13]. 
In nonlinear dissipative systems, localized modes appear 
due to an additional balance of gain and loss at some special 
values of parameters and the existence of families of local¬ 
ized modes is very rare [14,15]. While the PT-symmetric 
systems can be viewed as a special class of dissipative sys¬ 
tems, they are quite unique, combining many features of 
conservative and dissipative systems, being able to support 
families of localized and periodic modes. 


In this paper, we review the studies of nonlinear phe¬ 
nomena in optical PT symmetric systems. We discuss both 
theoretical concepts and predictions, as well as a relatively 
small yet quickly growing number of experimental demon¬ 
strations. We start with the discussion of the simplest con¬ 
figuration of PT-symmetric systems in optics in the form of 
a pair of coupled waveguides, an optical coupler, in Sec. 2. 
Then, in Sec. 3 we extend the discussion to multiple cou¬ 
pled waveguides, or oligomers. Further in Sec. 4 we analyse 
the properties of solitons in couplers and extended periodic 
systems. In Sec. 5, scattering of waves incident on nonlinear 
PT structures is considered. In Sec. 6, we discuss effects in 
PT structures, where gain and loss are modulated along the 
propagation direction. Finally, in Sec. 7 we present conclu¬ 
sion and outlook. 


2. PT-symmetric couplers 

A pair of coupled waveguides with gain and loss represents 
the simplest configuration of PT-symmetric optical system. 
However, such a system can already showcase rich physics 
and phenomena associated with the interplay of nonlinearity 
and PT-symmetry. In the following, we first overview the 
properties of linear couplers, and then discuss the effect 
of nonlinearity on the optical modes and the phenomenon 
of nonlinearly-induced PT symmetry breaking, as well as 
some generalizations. 


2.1. Linear properties 


Directional couplers composed of waveguides with gain and 
loss regions can be used to realize PT-symmetric optical 
structures. In the pioneering papers [4,5], different designs 
of optical structures were suggested, where optical beams 
would exhibit an effective PT-symmetric potential. In such 
structures, the real part of the optical refractive index should 
be symmetric, and the imaginary part - antisymmetric, i.e. 

n(x,y) =n*(-x,y). (2) 

This means that the absolute amounts of gain and loss should 
be the same, and gain/loss regions should have mirror con¬ 
figurations with respect to the central symmetry point. 

For experimental realization of optical PT couplers, di¬ 
rectional couplers with different values of gain or loss in 
each of the waveguides were considered [6,7], see Fig. 1(a). 
Beam propagation in such structures can be described by 
a pair of coupled equations for the amplitudes of modes 
in each of the waveguides, which can be represented in a 
Hamiltonian form as 



Ha, a(z) 


ai(z)\ H= f iP 1, ~C \ 
«2(z)/’ V ~C, -ipi) 


(3) 


where z is the propagation distance, a\ and a 2 are the mode 
amplitudes in the first and second waveguide, respectively, 
P 12 define the rates of gain in the first waveguide and loss 
in the second waveguide, and C is the coupling coefficient 
between the modes of two waveguides. 

The Hamiltonian possesses PT symmetry when applied 
together with Gauge transformation [7], 

PT(H - pi) = (H - pI)PT, (4) 

where parity operator swaps the mode amplitudes between 
the waveguides, P = (0,1; 1,0), T is a time-reversal operator 
which changes z — z and performs a complex conjugation, 
I is an identity matrix, and p = (pi — P 2 )/2 defines the 
average gain or loss between the two waveguides. 

Important features of PT couplers become apparent by 
studying the properties of optical supermodes. These so¬ 
lutions have the form a ^ m \z) = a^(0)exp(/j8 m z), where 
m = 1,2 is the mode index. The mode propagation con¬ 
stants are f3 m = Po + jS m , j6 m = Ccos(<p±) and the pro¬ 
files are a^ = exp(— i(p±)af l \ where sin(<p±) = p/C, 
P = (Pi +p 2 )A Po = i(Pi -P 2 )A 
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Figure 1 Linear PT couplers, (a) Characteristic real (n R , red line) and imaginary (nj, green line) parts of the complex refractive-index 
distribution, (b) Imaginary eigenvalue part for optical supermodes characterizing their gain or loss vs. the gain/loss coefficients in 
individual waveguides. (c,e) Calculated optical supermodes and dynamics below the critical gain/loss value: (c) the supermode intensity 
is equally divided between the two waveguides and (e) periodic beating of light between two waveguides, with zero average gain or 
loss, when the system is excited at either channel 1 or channel 2. (d,f) Regime above the critical gain/loss value: (d) a supermode is 
isolated at one waveguide and (f) light localizes in waveguide with gain and exhibits amplification. Figures (a,e,f) after Ref. [6]; (b,c,d) 
after Ref. [7]. 


When \p/C\ < 1, i.e. when the amount of gain and 
loss is below a critical value , then cp± are real, accord¬ 
ingly Irn/3 m = 0, and both supermodes exhibit the same 
amount of gain/loss as Im/3i = Im/32, see Fig. 1(b). In 
this regime the supermode profiles are PT-symmetric, i.e. 

PTa( m )(0) = exp [i(p± — 2/arg(«[”^)]a^(0). Accordingly, 
the mode intensity is symmetrically distributed between the 
waveguides, i.e. \a\ \ = \a2\, see Fig. 1(c). 

However, when the gain/loss exceeds a critical value, 
\p/c\ > 1, then (p m become complex, supermodes ex¬ 
hibit different gain/loss [Fig. 1(b)], and their profiles be¬ 
come asymmetric and no longer conform to PT symmetry 
[Fig. 1(d)]. Accordingly, this regime is called broken PT- 
symmetry, referring to the symmetry of the supermodes. 

An input optical beam in general excites a superposi¬ 
tion of the supermodes, which properties thereby determine 
the beam dynamics. Below the threshold, the beating of 
two supermodes can lead to a periodic switching of light 
between the coupled waveguides, see Fig. 1(e). Such be¬ 
havior is similar to a conventional conservative coupler, 
and indeed the power is conserved on average after every 
oscillation between the lossy and amplifying waveguides. 
The beating period is z p = 27r/|/3+ — j8_| = 7r(C 2 — p 2 ) -1 / 2 , 
and it grows monotonously as the gain/loss is increased, ap¬ 
proaching infinity close to the critical threshold. Above the 
threshold, one supermode decays, and the remaining su¬ 
permode primarily localized in the waveguide with gain 
determines amplification irrespective of the input condition, 
see Fig. 1(f). 


2.2. Nonlinear coupler and symmetry breaking 


The effect of nonlinearity on the beam dynamics in direc¬ 
tional couplers composed of waveguides with gain and loss 
was originally described theoretically already two decades 
ago in Refs. [11, 12]. It was predicted that such structures 
can offer benefits for all-optical switching in the nonlinear 
regime, lowering the switching power and attaining sharper 
switching transition. In the last decade, the renewed interest 
in such structures as realizations of symmetric optical 
systems has led to a series of extensive theoretical and more 
recently experimental studies. 

We first consider the effect of Kerr-type nonlinearity, 
which can modify the refractive index in each of the wave¬ 
guides, depending on the optical intensity [11, 12,16-19]. 
This effect can be modeled by including nonlinear terms in 
the coupled-mode Eq. (3) as follows [17], 


,da\ 

dz 

,da2 

dz 


r\ 

ipd\ Cci 2 + G(\a\ | )$i = 0, 

r\ 

— ipa2 + Ca\ + G(\a2\)a2 = 0 , 


(5) 


where the function G characterizes the nonlinear response 
(we assume it to be real-valued), and p > 0 is the loss/gain 
coefficient in the first and second waveguides. We will 
consider here the regime below the linear £?^-symmetry 
breaking threshold, when p < C. 

It was found that nonlinear solutions belong to two 
classes [17]: (i) periodic solutions, where the intensities and 
relative phases in two waveguides are exactly restored after 
each period (z z~\~z p ) [Fig. 2(a)], or (ii) solutions where 
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the total intensity grows without bound due to nonlinearly- 
induced symmetry breaking [Fig. 2(b)]. This classification 
is valid for arbitrary Kerr-type nonlinearities with smooth 
response functions G(7). Quite remarkably, in case of cubic 
nonlinear response with G(7) = y/, Eqs. (5) were shown to 
be integrable and their solutions can be formulated analyt¬ 
ically [16,18-22]. The dependence of the minimal input 
intensity (7 cr ) required for nonlinear switching on gain/loss 
coefficient is shown in Fig. 2(c), which illustrates that the 
threshold for nonlinear switching is drastically reduced for 
larger gain/loss coefficients. 

It is convenient to represent the mode amplitudes in the 
following form, 

Cl\ = V / /(z)cos[0(z)]exp[+i'(p(z)/2] exp[i/}(z)], 

( 6 ) 

«2 = y7(z)siii[0(z)]exp[-/<p(z)/2]exp[/j3(z)], 

where I is the total intensity, 0 and (p define the relative 
intensities and phases between the two input waveguides, 
and jS is the overall phase. Then, the initial conditions corre¬ 
sponding to different solution types (oscillating or growing) 
can be conveniently visualized on a phase plane, see an 
example in Fig. 2(d). This plot illustrates that even for high 
intensity above the threshold, the type of dynamics depends 
on the relative amplitudes and phases in the two waveguides. 
Interestingly, the type of nonlinear dynamics remains the 
same if we swap the intensities between the two waveguides, 
which corresponds to a transformation G —^ tt/2 — 0. In par¬ 
ticular, we could couple light at the input just to the wave¬ 
guide with loss, or to the waveguide with gain, and the type 
of dynamics would be the same. This is a counter-intuitive 
result, since in the first case the total intensity will initially 
decrease, whereas in the second case the total intensity will 
be growing. However, in both cases the type of dynamics 
will be determined only by the initial intensity level. This is 
a highly nontrivial consequence of linear PT-symmetry in 
the strongly nonlinear regime. 

Whereas in general solutions are oscillating or grow¬ 
ing, for particular initial conditions wave intensities can 
remain constant. Such regime corresponds to the excita¬ 
tion of stationary nonlinear modes, which can be viewed 
as a generalization of linear supermodes. They are posi¬ 
tioned at the phase space points (7 = 7o,0 = ;r/4,<p = <p±) 
and j6 = (5±z , where f3± = G(7o/2) +Ccos(<p±), sin(<p±) = 
p/C, and cos(<p±) = =p[1 — (p/C) 2 ] 1 / 2 . The fixed point in 
phase space is a stable center if (i) y < y cr and (p = <p_ or 
(ii) y > — y cr and (p = <p + , where y cr = | cos(<p±)| and y = 
G / (7o/2)7o/(2C) (prime stands for the derivative). If these 
conditions are not satisfied, then the fixed point is a saddle, 
which indicates an instability. In case of self-focusing Kerr 
nonlinearity, G(7) = y/ with y > 0, the stationary point at 
<p + is always stable at arbitrarily high intensities, whereas at 
<p_ the instability appears for 7 q > 7o C r = 2[1 — p 2 ] 1 / 2 . The 
stable point at <p + is indicated in Fig. 2(d) by a star, and un¬ 
stable at <p_ - by a triangle. We note that the unstable point 
is located at a boundary between the domains of periodic or 
growing solutions. 



P q> 


Figure 2 Nonlinear PT coupler. (a,b) Intensity dependencies on 
propagation distance in the first (dotted line) and second (dashed 
line) waveguides; the solid line shows the sum of individual intensi¬ 
ties. Shown are the regimes of (a) periodic oscillations and (b) non¬ 
linear localization and amplification in the waveguide with gain, 
(c) Minimum critical intensity required for nonlinear PT-symmetry 
breaking vs the gain/loss coefficient, (d) Regions of PT symmetry 
(white shading) and symmetry breaking with nonlinear switching 
(black shading) in the plane of initial conditions. Stars and trian¬ 
gles mark stable and unstable stationary solutions, respectively. 
Parameters are G(7) = y/, y = 1, C = 1, (a,b,d) p = 0.5 and 
I = 2.2, (a) (p = 7i/6 — 7t/20, (b) (p = tz/6 + 7t/20. Figures after 
Ref. [17]. 

Nonlinear PT couplers can have applications for effi¬ 
cient all-optical signal manipulation with reduced switching 
power and the possibility to perform intensity-dependent 
amplification. Additionally, such couplers can operate as 
unidirectional optical valves [16] when nonlinearity breaks 
PT symmetry, as in this case the output beam becomes 
localized in the gain channel, irrespective of initial condi¬ 
tions. It was suggested that switching performance can be 
enhanced with application of Bragg gratings, involving spe¬ 
cially engineered modulation of the real and imaginary part 
of the optical refractive index [23]. These concepts can be 
extended to other physical systems beyond optics, including 
electrical circuits [24,25]. 

A PT-symmetric coupler, with additional gain and 
loss proportional to nonlinear terms has been studied in 
Refs. [26-28]. The system can be described by the follow¬ 
ing equations, 

r\ 

iu\ = -U 2 + (x-iYnon)\u\\ u\-iyu\, (7) 

r\ 

iu 2 = -Ml +{x + iynon)\u2\ «2 + *y«2- (8) 

where % determines the nonlinearity strength and y non ac¬ 
counts for the PT-symmetric nonlinear loss and gain. The 
symmetric and asymmetric eigenstates of the system are 
found analytically. The additional nonlinear gain and loss 
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result in the emergence of asymmetric solutions which not 
only are involved in symmetry-breaking bifurcations, but 
also produce asymmetric linearization matrices with spectral 
properties that reflect this asymmetry. If such nonlinear PT- 
symmetric coupler implemented into a linear chain, it may 
give raise to new types of nonlinear Fano resonances, with 
entirely suppressed or greatly amplified transmission [27] . 


2.3. Actively coupled waveguides 


Nonlinearity can provide an important mechanism to bal¬ 
ance gain and loss on average [22,29,30]. Such situation 
can arise even when linear PT symmetry is absent, in par¬ 
ticular when the coupling between the waveguides is non¬ 
conservative and provides gain or loss depending on the 
relative phase between the guided modes [29] . Physically, 
this can be realized when a pair of nonlinear optical wave¬ 
guides with absorption are placed in a medium with power 
gain. The evolution of mode amplitudes in two waveguides 
is governed by the following normalized equations, 

d Vf\ i 1 9 / x 

i—- + iyy i +1 v'i I Vi + ( 1 - ia ) V2 = o, 

* (9) 

dxif 2 I ,9 /.X 

i—— +iW2 + \W2\ W2 + (l ~ia)Vi =0, 

dz 

where y > 0 is the net loss rate and a is the active coupling 
coefficient. Note that these equations are different from 
Eqs. (5), and it was found that nonlinear dynamics demon¬ 
strates distinct regimes including stationary, periodic, and 
chaotic oscillations as outlined in Fig. 3(a). 

The reason for the appearance of oscillations is the fol¬ 
lowing. For a > y, symmetric supermodes with i ff\ — \f /2 
exhibit a linear amplification rate (y— a). As the amplitude 
increased beyond a certain level, nonlinearity results in sym¬ 
metry breaking and localization of light in one waveguides, 
and then the mode exhibits linear loss with rate approaching 
y. Therefore, real nonlinearity acts as an effective nonlinear 
absorber. Examples of the resulting oscillations are shown in 
Figs. 3(b,c). An important advantage of the actively coupled 
nonlinear PT structure is its structural stability. For the given 
loss rate, the system supports stationary and periodic light 
beams in a wide range of gain coefficients. This is a fun¬ 
damental distinction from the conventional PT-symmetric 
coupler discussed above, where one has to tune the gain to 
match the loss exactly to obtain periodic oscillations. 

Another type of PT symmetry can appear when active 
coupling between the modes is nonlinear [31]. This can 
be potentially realized in a two waveguide coupler in the 
presence of a nonlocal nonlinearity, which respects PT sym¬ 
metry. The corresponding model shows unique behavior, in 
particular admitting both bright and dark soliton solutions 
at the same time. 



Figure 3 Actively coupled PT system: two lossy waveguides 
placed in an active medium with gain, (a) The chart of attractors 
in the plane of the normalized net loss rate (y) and active coupling 
coefficient (a), (b) The limit cycle with two oscillations per period 
(a = 9, y= 4.2). (c) The chaotic attractor (a = 9, y = 7.6). Figures 
after Ref. [29]. 

2.4. Anti-PT-Symmetry and parametric 
amplification 

Nonlinear modes in PT potentials can also be supported 
by quadratic (x^) optical nonlinearity [32, 33]. Com¬ 
pared to the case of Kerr-type nonlinearity discussed above, 
quadratic nonlinear response involves parametric coupling 
of the fundamental and second-harmonic optical waves. 
Such interactions give rise to a rich family of localized 
solitons. 

Quadratic nonlinearity can also facilitate a fundamen¬ 
tally important regime of parametric amplification, which 
is efficiently realized through difference-frequency genera¬ 
tion [34]. Here the amplification rate is determined by the 
pump, enabling ultrafast all-optical tunability. Quadratic 
nonlinear PT-symmetric systems can, on one hand, realize 
ultrafast spatial signal switching through pump-controlled 
breaking of PT symmetry, and on the other hand enable 
spectrally-selective mode amplification [35]. The process 
of optical parametric amplification in PT coupler based 
on nonlinear mixing between a strong pump, and signal 
and idler waves is illustrated in Fig. 4(a). The evolution 
of near-degenerate (with close similar frequencies) signal 
and idler waves in the undepleted pump regime is modelled 
by coupled-mode equations, which can be represented in a 
Hamiltonian form: 

5a 

i— = Ha, (10) 

dz 
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signal—^ 
pump 
signal—^ 



(a) 




1 




Figure 4 PT parametric amplifier: (a) scheme of wave mixing 
in a quadratic nonlinear coupler with linear absorption in one 
waveguide, (b) Number of spectrally anti-PT symmetric mode 
pairs, (c) the largest mode gain (white line marks zero), (d) fraction 
of signal intensity in the first waveguide at the output vs. the input 
pump amplitude in the first waveguide and the phase-mismatch. 
Figures after Ref. [35]. 


T 

where a(z) = (a s i(z),a s 2 (z),a* l (z),a* 2 (z)) is a vector of 
wave amplitudes in the two waveguides, the subscripts stand 
for signal (‘s’) and idler (T) waves in two waveguides (‘1’ 
and ‘2’), 


fP-iYi ~C iAi 0 \ 

H= -C p - iyi o iA 2 (11) 

iA\ 0 -P-iyi C 

\ 0 iA* 2 C -P-iyiJ 

jS defines the phase mismatch of parametric wave mixing, y 
are the linear loss coefficients, C is the mode coupling coef¬ 
ficient between the waveguides, and A are the normalized 
input pump amplitudes. 

Remarkably, the Hamiltonian possesses a spectral anti- 
PT symmetry , corresponding to a negative sign on the right- 
hand side of the following equality, 


Characteristic dependance of the number of PT-symmetric 
mode pairs in depending on the pump amplitude in the first 
waveguide and the phase-mismatch is shown in Fig. 4(b), 
and the largest mode gain is presented in Fig. 4(c). Re¬ 
markably, the established relations of mode symmetry and 
gain/loss are reversed in comparison to previously discussed 
spatial PT-symmetry in directional couplers, due to the spec¬ 
tral anti-PT symmetry of parametric wave mixing. Specif¬ 
ically, when all modes have broken spectral PT symmetry 
[blue shaded regions in Fig. 4(b)], the pairs of eigenmodes 
exhibit the same amount of gain/loss, and effectively the 
amounts of gain and loss are averaged out between the 
eigenmodes. However, upon transition to the region with 
spectrally PT-symmetric modes [green and red shaded re¬ 
gions in Fig. 4(b)], there appears an unequal redistribution 
of gain and loss between the modes. One PT-symmetric 
eigenmode exhibits gain much larger then all other modes, 
while the latter experience stronger loss. Such sensitivity 
of amplification to PT-breaking threshold could be used to 
discriminate between multiple spectral modes, analogous to 
the concept of PT-lasers [8,9]. 

In case of perfect phase-matching, jS = 0, and pump am¬ 
plitudes of the same phase Im(Ai) = Im(A 2 ) = 0, the Hamil¬ 
tonian also features spatial PT symmetry. Both the spatial 
and spectral PT-symmetry breaking occurs simultaneously 
at the threshold |7i — 72 — (Ai + A 2 ) | = 2C. However, the 
spatial and spectral symmetries are opposite: a mode pair is 
spatially PT-symmetric and has spectrally broken symmetry 
below threshold, whereas the situation is reversed above the 
threshold. Numerical simulations reveal a strong connection 
between spectral symmetry and spatial dynamics even for 
non-zero phase mismatches. Indeed, an increase of pump 
amplitude can control the period of mode coupling between 
the waveguides, while the oscillations get suppressed close 
to the spectral PT threshold, see Fig. 4(d). It is expected that 
due to the universality of parametric amplification processes, 
these concepts can be extended to different physical mecha¬ 
nisms including four wave mixing in media with Kerr-type 
optical nonlinearity. 


2.5. Nonlinear modes 


P 1 ? +P 2 ,+TH = —HPi + P 2 + T. (12) 

Here T is a usual time-reversal operator which changes 
z —^ — z and performs a complex conjugation. The parity 
operators operate in spectral domain , interchanging the 
signal and idler waves, 

Pi,± = {a s \ ^ ± 0 * 1 } ? P 2 ,± = {a S 2 ^ ± a a} • (13) 

The Hamiltonian is linear, and the properties of its solutions 
are defined by the spectrum of its four eigenmodes. Over¬ 
all, there can be three possible symmetry regimes: (i) two 
mode pairs with broken spectral anti-PT symmetry, (ii) one 
pair of PT-broken modes and a pair of anti-PT symmet¬ 
ric modes, or (iii) two pairs of anti-PT symmetric modes. 


Analysis in previous sections was based on the consideration 
of coupling between the fundamental guided modes. How¬ 
ever, new effects can appear when PT systems support mul¬ 
tiple modes, including bound and extended (non-localized) 
states. Various types of double-well potentials have been 
investigated [36-38]. It was confirmed that discrete cou¬ 
pled equations can accurately describe the full continuous 
dynamics in a broad range of parameters. However, contin¬ 
uous models can be used to predict a range of important 
effects beyond two-mode approximation, such as the appear¬ 
ance of radiation losses in curved PT waveguides [36]. On 
the other hand, it was predicted that a PT structure placed 
in a nonlinear medium can trap soliton [39]. Furthermore, 
soliton scattering by a localized mode of PT coupler can ex¬ 
hibit nonreciprocal and non-conservative scattering, which 
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is strongly dependent on the relative phase between the 
mode and the soliton [40]. 

Properties of nonlinear modes supported by a harmonic 
PT-symmetric potential also show a number of distinguish¬ 
ing features [41]. In particular, the modes bifurcating from 
different eigenstates of the underlying linear problem can 
belong to the same family of nonlinear modes. By adjusting 
the potential profile, it is possible to enhance the stability of 
small-amplitude and strongly nonlinear modes compared to 
the case of conservative real harmonic potential. 

3. Discrete PT-symmetric oligomers 

PT-symmetric discrete oligomers have recently attracted a 
lot of attention, motivated by possibilities of their exper¬ 
imental realisation [ 6 , 16]. In particular, many scientific 
works address the questions of existence and stability of 
nonlinear states in PT-symmetric oligomers, which may 
drastically differ from eigenmodes of underlying linear sys¬ 
tems. The nonlinear effects in PT-symmetric systems can 
be utilized for an efficient control of light including all- 
optical low-threshold switching and unidirectional invisi¬ 
bility [16,42]. The possibility to engineer PT-symmetric 
oligomers, which may include nonlinearity, triggers a broad 
variety of studies on both the few-site systems and en¬ 
tire PT-symmetric lattices, including one-dimensional PT- 
symmetric dimer [27,43], trimer [43,44], quadrimer [43,45], 
2D PT-symmetric plaquettes [45,46], and PT-symmetric fi¬ 
nite/infinite chains [47,48] and necklaces [49] and multicore 
fibers [50]. 

In this section we discuss PT-symmetric configurations 
consisting of three or more coupled waveguides. We first 
formulate the general model of finite discrete PT-symmetric 
oligomer and next consider some special configurations of 
PT-symmetric trimers and quadrimers. We also discuss finite 
open and periodic few-site geometries. 

Consider an array of N waveguides coupled to each 
other. We assume that each waveguide experiences some 
gain or loss determined by parameter y n , where n is a wave¬ 
guide number. Taking into account conservative Kerr nonlin¬ 
ear response, the light propagation through this system can 
be modeled by a system of discrete nonlinear Schrodinger 
equations (DNLSE), 

N _ 

illn — ^ KnmUm \Wn\ M n ^Yn^n> M — 1 • N , (14) 

m— 1 

where u n are mode amplitudes, u n = du n /dz, z is the propa¬ 
gation distance, y n > 0 (y n < 0 ) corresponds to loss (gain), 
and K nm are the coupling coefficients between sites n and 
m. We assume that K = {K nm } is a real symmetric matrix, 
i.e. K nm = K mn = K* m . The system (14) can be rewritten in 
matrix form as 

m = — [H + F(u)]u, H = K-WG, (15) 

where G = diag(yi, y 2 ,--,Yn) 

and F(u) = diag(|wi | 2 , | w 2 1 2 ? •••, | | 2 ) • (16) 



Figure 5 Schematic of PT-symmetric trimer for (a) open and (b) 
periodic boundary conditions. Signs and ”+” denote loss and 
gain sites, respectively, k - coupling coefficient. 

Note that if H commutes with the PT operator, i.e. 
[PT,H] = 0, then H is a PT-symmetric matrix. 

Stationary nonlinear modes of Eq. (15) are sought in the 
form u(z) = exp(/fe)w, with propagation constant (eigen¬ 
value) b and w satisfying the following equation 

bw = [H + F(w)]w. (17) 

N 

Hereafter we assume that £ y n = 0, which is a necessary 

n —1 

condition for the system spectrum to be real [45]. 

In the following, we consider in detail two PT-symmetric 
geometries with the nearest-neighbor interaction and ho¬ 
mogeneous coupling strength between the sites. Then the 
coupling coefficients can be expressed as K nm = kS\ n _ m \ \. 


3.1. PT-symmetric trimer 


PT-symmetric trimer, first discussed in [43] and revisited 
in [26,51] is schematically shown in Fig. 5. It is com¬ 
posed of three waveguides: active, lossy and conservative 
ones. In the linear regime there is only one possible PT- 
symmetric interposition of these waveguides with non-zero 
PT-symmetry breaking threshold: equal gain and loss pa¬ 
rameters, i.e. Y\ = — 73 = y, with conservative waveguide in 
the middle. The PT-symmetry breaking threshold is defined 
by the relation y cr j t = y/2k. 

In the nonlinear regime, the dynamical equations for the 
open trimer [open boundary conditions, Fig. 5(a)] are 

i 19 

iu\ = — ku 2 ~ \ui\ u\—iyu\, 

r\ 

iU2 = — k(u\ + M 3 ) — \U2\ U2, (18) 

1 19 

iu?, = — ku 2 — | M 31 M 3 + iyu 2 >. 

Here we are interested in the interplay of nonlinearity with 
PT-symmetry and its effect on stationary modes and their 
stability. We search for the stationary solutions of Eqs. (18) 
in the form u\ = Aexp[ibz + <f> a ], U 2 = B exp [ibz + <fe], and 
M 3 = Cexp[/fe + (/) c ]. Substituting these ansatz into Eqs. (18) 
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we find relations on the soliton parameters, 

A 2 [y 2 + (/? — A 2 ) 2 ] 2 — k 2 b [y 2 + (b—A 2 ) 2 - 


2 k 4 A 2 + 2 k 4 b = 0, 

2 Z>±^Z> 2 -8A 2 (Z>-A 2 ) 

B - 


C = A. 


sin(<fe - <p a ) 

COS ((/)/? - <j ) a ) 


- sm(<p b - <p c ) 

cos (<j) b - (pc) = 


yA 

M-A 3 

kB 


(19) 


Note that without loss of generality we may put = 0, 
which yields A = C,(j) a = —(j) c and thus 


u 3 



Thereby, all stationary modes of the PT-symmetric trimer 
possess equally distributed intensity between loss and gain 
waveguides. 

Analysis of Eqs. (19) shows that there are up to five (at 
least one) branches of stationary solutions of Eqs. (18) de¬ 
pending on the particular values of b and k. The dependence 
of these branches on gain/loss parameter y for particular 
case b = k = 1 has been discussed in detail in [43, 51]. 
Three different branches were found: first branch is mostly 
unstable (except small region that splits unstable region into 
two domains), while the second one is chiefly stable. These 
branches collide in a saddle-center bifurcation at y = 1.043 
and disappear. Interestingly, the third branch bifurcates from 
zero for y > y /2 — b 2 and persists for larger values of y even 
beyond the PT-symmetry breaking threshold y cr i t = y/2, al¬ 
though this branch is unstable for y > 1.13. In Fig. 6 the 
system dynamics for the initial conditions belonging to the 
first branch is presented for two different values of gain/loss 
parameter y. In Fig. 6(a) for y = 0.5 (first unstable region), 
the nonlinear modes are unstable and the growth of intensity 
is observed in the active site u 2 , while two other modes u\ 
and u 2 decay. In Fig. 6(b) for y = 1.04 is taken from the 
second unstable region, but this time intensity growth is 
observed not only in the active waveguide but also in the 
conservative one. Both these types of dynamics are also 
observed for the second branch. 

Note that in the considered case it is possible to find two 
additional branches of solutions of Eqs. (19) assuming b to 
be complex-valued, although modes with such propagation 
constants will not be longer solutions of the original sys¬ 
tem (18). Such so-called ghost states can predict evolution 
of the systems (18) at short propagation distances [51]. 

The second configuration of the PT-symmetric trimer 
schematically depicted in Fig. 5(b) has also been investi¬ 
gated in Ref. [43] (and also discussed in [52]). Although 
the PT-symmetry breaking threshold for underlying linear 
system is y cr j t = 0, the nonlinear branches of solutions still 
exist for nonzero values of gain/loss parameter y. Thus, 
nonlinearity provides stable stationary nontrivial solutions 
of the system, the linear spectrum of which has complex 
eigenvalues. 



Figure 6 Dynamical evolution of initial data belonging to the 
first branch of the PT-symmetric trimer. (a) y = 0.5 (first unstable 
region), (b) y= 1.04 (second unstable region). Red, blue and 
green lines denote w 3 , u 2 , and zq, respectively. Adopted from [43]. 
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Figure 7 Schematic of PT-symmetric quadrimer with nearest- 
neighbor interaction with equal coupling strength. 


More complicated structure of PT-symmetric trimer in¬ 
cluding competing nonlinear gain with linear loss in the first 
waveguide and vise versa in the third waveguide has also 
been studied in [26]. The introduced nonlinear gain and loss 
give rise to rich features of asymmetric modes not observed 
in linear gain/loss profile systems. 

It is worth to note here that experimental realization of 
dissipative trimer having ’’gain-loss-gain” profile has been 
done in [51]. Although this trimer does not possess PT- 
symmetry, this experimental setup can be considered as a 
first step to realization of PT-symmetric trimer. 


3.2. PT-symmetric quadrumer 

In contrast to the PT-symmetric trimers considered above 
the four-site geometry admits several PT-symmetric con¬ 
figuration. Taking into account PT-symmetry requirements 
(symmetrical distribution of gain and loss sites) the general 
model of a one-dimensional PT-symmetric quadrimer with 
nearest-neighbor interaction [see Fig. 7] can be described 
by the following equations, 

1 1 9 

iii\ = — ku 2 — \u\\ u\ — iy\u \, 

r\ 

iu 2 = —k(u\ + m 3 ) - \u 2 \ u 2 -iy 2 u 2 , ( 20 ) 

r\ 

iu 2 = — k[u 2 + W 4 ) — \u 2 \ u 2 -\-iy 2 U 3 , 

1 1 9 

iii\ — —ku 2 — \U 4 \ U 4 + ij\U 4 . 


Similar to the trimer case, the existence and stability of 
stationary nonlinear modes for yi = Ji has been investigated 
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b b b 


W 2 = W3 which leads to the set of equations 

bw 1 =W2 + |wi| w\ + ij\W \, (22) 

bw2 = w\ +W2 + |w2| 2 W2 + /72W2. (23) 

Substituting vtq = Wiexp(/<I>) and W2 = W2exp(/d>) with 
real W\ and into Eq. (22), we get W2 = W\(b — Wf — 
iy ), while from Eq. (23) we have exp(—2/<f>) = f{W\) = 
(MT 2 - Wi - |W 2 | 2 W2 - iyW 2 )/W :2 *. The roots of the equa¬ 
tion \f(W \)| 2 = 1 immediately give us wi and W2. The latter 
equation can be reduced to the P$(Wf) = 0, where P% is a 
eighth-degree polynomial with real coefficients, each real 
root of which corresponds to a nonlinear branch of modes. 
Since the roots depend on b continuously, we have continu¬ 
ous families of nonlinear modes for fixed system parameters. 
Note that these families of modes are additional to one found 
in [43] for the case y\ =72. The branches of modes can be 
characterized by the total energy flow 


Figure 8 Regions of stability underlying linear system of 
Eqs. (20). The dark diamond-shaped region corresponds to unbro¬ 
ken PT-symmetry phase; the light-grey regions represent parame¬ 
ter ranges where two of eigenvalues are real and two are complex; 
and the white regions correspond to all eigenvalues are complex. 
Curves in (a)-(e) correspond to the branches of nonlinear modes, 
where U is a total energy flow determined by (24). The solid blue 
and dashed red lines indicate the stable and unstable branches 
of modes, respectively. Panels (a)-(e) correspond to the points 
(a)-(e) in the Phase Diagram. Adopted from Ref. [45]. 


in [43]. Such structure additionally possesses completely 
asymmetric solutions and the symmetry-breaking bifurca¬ 
tion structure is more complicated. 

Let us consider in detail general case of bi-parameter 
gain and loss, i.e. 71 7^ 72 [45]. The analysis of the underly¬ 
ing linear system of (20) shows that there are three types of 
regions in the parameter plane 7172 [see Fig 8 left top panel]: 
(i) unbroken PT-symmetry, with all eigenvalues real; (ii) 
broken PT symmetry with two real and two complex conju¬ 
gated eigenvalues (this case is only possible for 71 7^ 72); (iii) 
broken PT-symmetry with all eigenvalues complex. Note 
that there are ’’triple” points in left top panel Fig. 8 denoted 
Tj, which correspond to contact points of three regions. At 
these points, all eigenvalues b n = 0, hereafter ’’tilde” indi¬ 
cates eigenstates of the underlying linear system. 

Next, we discuss an effect of nonlinearity on stationary 
modes, considering unbroken PT phase where eigenstates 
w of the matrix H of the system (20) are simultaneously the 
eigenstates of PT operator, i.e PoTw = w up to an irrelevant 
phase shift. Here Pq is parity operator of the form 


^0 = 




( 21 ) 


Hereafter 01,2,3 are Pauli matrices and 0 is 2 x 2 zero matrix. 
It is reasonable to look for nonlinear modes possessing the 
same property = w. Thus, we may put w\ = w4 and 


1 4 

U(b) = 2 £ Kl 2 . (24) 

n =1 

In Fig. 8(a)-(e) the branches of solutions are depicted for 
the gain/loss parameters 71 and 72 corresponding to the 
points (a)-(e) in Fig. 8 left top panel. It is found that: (i) 
for small values of gain/loss parameters 7^2 from the exact 
PT-symmetry phase region, nonlinear modes bifurcate from 
linear ones; (ii) if 71.2 are chosen from an unbroken PT- 
symmetry region, but close to triple points 7) as shown in 
Figs. 8 (b), (c), nonlinear modes bifurcate from linear limit, 
but rapidly loose stability and two of these branches cease 
to exist as U reaches some critical value; (iii) the stable 
nonlinear modes can be found even in the region of broken 
PT symmetry [see Figs. 8 (d), (e)]. 

We now give a general perspective on the role of nonlin¬ 
earity for PT-symmetric systems [45]. It is known that linear 
PT-symmetric systems exhibit similar behaviour (e.g. real 
spectrum and energy conservation of eigenmodes) as conser¬ 
vative ones, if they operate below critical point. Thus, there 
exists a unitary equivalence between PT-symmetric and Her- 
mitian systems [53,54]. In terms of considered oligomers 
(we return to matrix representation (15)), it means that there 
is a unitary matrix R that transforms exactly PT-symmetric 
matrix (in PT unbroken phase) Hq into Hermitian one with 
the same set of eigenvalues, i.e. RHqR -1 = H# = H^. This 
new matrix H h will include coefficients y n , which will de¬ 
termine the coupling strength. Thus, the following question 
naturally arises: are there any significant differences be¬ 
tween these PT-symmetric and Hermitian systems having 
the same eigenvalues in nonlinear regime? 

To answer this question we introduce new system, de¬ 
scribed by the equation m = — [H# + F(u)]u, where H h 
can be constructed explicitly, following [53,54]. Families 
of stationary nonlinear modes of this new system were ob¬ 
tained in [45] and are presented in Fig. 9 for the same val¬ 
ues of 71,2 as in Figs. 7 (a)-(c). While for 71 2 taken from 
the ’’middle” of the unbroken PT phase the behaviour of 
mode branches for PT-symmetric and Hermitian systems 
looks similar (except their stability regions)[cf. Fig. 8(a) and 
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Figure 9 Families of nonlinear modes of the Hermitian quadrimer 
determined by H H . Parameters 7 ^ 2 , which now affect on the cou¬ 
pling strength, are the same as in panels (a-c) of Fig. 8. Adopted 
from [45]. 


Fig. 9(a)], for 772 placed near the border of the unbroken PT 
phase, the difference becomes significant [cf. Figs. 8 (b),(c) 
and Figs. 9 (b),(c)]. While for the Hermitian system the 
mode branches do not change qualitatively depending on 
772 and there is always at least one stable branch, in the 
PT-symmetric case there is a bifurcation of two branches 
and loss of stability for even small values of the energy flow 
U. 

From the result above, there could be an impression that 
an existence of continuous families of the nonlinear modes 
is an attribute of PT-symmetric systems. To check this con¬ 
jecture, we consider the case 71 = 72 = 7, but this time 
take into account not only the nearest-neighbor interactions, 
according to the following matrix: 


0 

1 

0 0 \ 

1 

0 

cos /3 — sin /3 

0 

cos/3 

sin 2/3 cos 2/3 

Vo 

— sin/3 

cos 2/3 — sin 2/3 / 


(25) 


where /3 determines a coupling strength. This structure is 
schematically shown in Fig. 10 in the left top corner. This 
configuration is symmetric with respect to the operator 


m) = (Jn P T)’ <26) 

where p (/ 3 ) = cos jS <J\ + sin /3 03. Note that Ho ( 0 ) coincides 
with the case considered above for 71 =72. 

It can be shown that Hq(/3) admits continuous families 
of nonlinear modes only for discrete number of points jS^ = 
nk/2 [45]. Thus, continuous families of nonlinear modes 
do not always exist for PT-symmetric system. 

Interestingly, by using appropriate transformation, e.g. 


M = 



(27) 


with /i = 03 + iO \, it is possible to obtain new PT-symmetric 
matrix Hi(/3) = MHo(/3)M -1 with the same set of eigen¬ 
values. Thus, the PT-symmetric quadrimer can be recon¬ 
figured to a PT-symmetric plaquette. In Fig. 10 the graph 
representations of both PT-symmetric quadrimer Ho(/3) and 
PT-symmetric plaquette Hi (jS) are shown for a generic case 
of arbitrary /3 and a particular one with /3 = %j 4. Here 
k± = \/2/2sin(jS d= tt/4). 

The PT-symmetric plaquettes can also be considered as 
blocks for construction of two-dimensional PT-symmetric 



Ho (/?) 

H x 09) 


1/1^ ^2 ^^3 U 4 

ztn 1 ^cos/^cos 2/3^ 

“1 n 

D~ k -( 

jh 

generic /3 

k+ 

w 

Y r 

k ._ 

tt7- 

— sin /3 

/3 = 7t/4 

i/i ^ i/i 2 

1 V 2/2 

m rh <0 (O 

u, 

( 

V 2/2 

u 4 C 

9 vj 

~k+ 

1 / 

u u 2 

u 3 

-V 2/2 

^-©» 

-V2/2 2 


Figure 10 PT-symmetric quadrimer H 0 (/3) and the correspond¬ 
ing plaquette configuration Hi(/3) obtained by means of M trans¬ 
formation of H 0 (/3). Adopted from [45]. 



Figure 11 (a) PT-symmetric birefringent arm coupler with prin¬ 

cipal optical axes of the two waveguides are tt/4 rotated with 
respect to each other, and (b) its graph representation in form of 
PT-symmetric plaquette. Red (”+”) and blue (”-”) colors denote 
active and lossy waveguides, respectively. Adopted from [55]. 


lattices. Therefore, it is important to figure out what are the 
properties of the nonlinear modes in such structures. Four 
types of fundamental PT-symmetric plaquettes correspond¬ 
ing to different interpositions of active and lossy waveguides 
were studied in Ref. [46] (one of possible realisation of such 
a plackets is presented in Fig. 1 1(b)). Existence of solution 
branches past the critical point of the underlying linear sys¬ 
tem and the appearance of additional asymmetric modes 
were reported, although most of the stationary modes were 
found to be unstable. 

The plaquette models can describe real physical sys¬ 
tems such as PT-symmetric multi-modes fibers [55, 56], 
including PT-symmetric coupler with birefringent arms was 
considered in [55]. The schematic of such a coupler and 
its graph representation are shown in Fig. 11. The coupling 
between elements includes both linear and nonlinear nature, 
which complicates the mode structure. Stationary modes 
with different polarisation were obtained and regions of their 
stability were identified. 


3.3. Multicore fibers 

Being driven by rapid development of multicore optical 
fibres, the theory of coherent propagation and power transfer 
in low-dimensional array of coupled nonlinear waveguides 
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Figure 12 Schematic (a),(b) radially symmetric multicore fibers 
with balanced gain and loss. From Ref. [50], and periodic PT- 
symmetric configurations: (c) N - site ring with embedded PT- 
symmetric coupler. Phase circulation direction of vortices with 
topological charge m > 0 and m < 0 is shown by anticlockwise 
and clockwise arrows, respectively; (d) PT-symmetric necklaces 
of alternating and clustered gain and loss waveguides. Adopted 
from [49,52], respectively. 


has been suggested [57,58]. This offers an opportunity to 
explore the multicore waveguide systems with PT-symmetry 
in the presence of gain and loss. Existence, stability, and 
dynamics of both linear and nonlinear stationary modes 
propagating in radially symmetric multicore waveguides 
with balanced gain and loss have recently been studied [50]. 
A multicore waveguide array of N identical waveguides 
arranged in a circular geometry with the central one being 
gain (loss) and all peripheral ones being loss (gain), as 
shown in Fig. 12(a)-(b), has been investigated. Taking into 
account only nearest-neighbor interaction, the evolution 
equations for the mode amplitudes read as, 

dA N 

—— = (fib + Oh) A + Co T B I + G|A| 2 A, (28) 

dz j=l 

dB • 

= (£i + ij\ )Bj +Ci (Bj+ 1 -\-Bj- 1) +CoA + G\Bj\ 2 B h 

where A, £q and yo denote the mode amplitude, refractive 
index, and gain/loss parameter in the central core waveguide, 
respectively, while Bj , £\ and are the same but in the j- th 
waveguide on the ring. Here Bq = B^ and B^+i=Bi and the 
nonlinear coefficients for all waveguides are identical and 
characterized by G. Co is the coupling coefficient between 
the core and other waveguides on the ring, while C\ stands 
for the coupling between two nearest-neighbors on the ring. 

It has been demonstrated that the system can generally 
be reduced to an effective PT-symmetric dimer with asym¬ 
metric coupling. In the linear regime, G — 0, the bounded 


dynamics can be observed in the limit of the effective dimer 
and for the active core waveguide. Here two modes with real 
propagation constants exist before a critical value of gain 
and loss for an onset of PT-symmetry breaking, while other 
modes always have imaginary propagation constants. This 
directly relates to a stable (unstable) propagation of modes 
when gain is localized in the core (ring) of the waveguide 
structure. 

In the case of self-focusing nonlinearity (G > 0), an 
interplay between nonlinearity, gain and loss results in a 
high degree of instability, and quasi-stable propagation may 
appear under some parameter conditions. By reducing to 
the effective dimer, several stationary modes were found, 
however these stationary modes turned out to be unstable 
in the presence of nonlinearity and nonzero gain and loss. 
More interestingly, when a spatially periodic modulation is 
applied, stabilization of nonlinear modes can be achieved. 

Nonzero energy flow between gain and loss waveguides 
may affect such entities as discrete vortices in the ring geom¬ 
etry including PT-symmetric elements [see Fig. 12(c)] [52]. 
One of the important characteristics of the vortex modes is 
a topological charge m which is equal to the phase winding 
number in units of 2k around the vortex origin and the sign 
of m determines the direction of power flow. It has been 
shown that degenerate linear vortex modes spontaneously 
break PT symmetry, however nonlinear propagation modes 
still can be found. By breaking T symmetry, the existence, 
stability and the dynamics of nonlinear vortex modes be¬ 
come sensitive to the sign of their charge m, offering an 
additional degree of freedom for all-optical control of dis¬ 
crete vortices. 

Generalization of the two-channel dispersive cou¬ 
pler [59] to a PT-symmetric arrangement of 2 N dispersive 
waveguides has been done in [49]. The arrays of wave¬ 
guides in the form of periodic necklaces with alternating 
and clustered gain and loss distribution [see Fig. 12(d)] 
were discussed and asymptotic behaviour of the PT break¬ 
ing threshold with number of sites N was derived. Thus, 
for the alternating array y cr j t ~ 1 /TV, while for the clustered 
geometry y cr i t ~ 1/A 2 as N «>. 

Note that, the questions of existence, stability and dy¬ 
namical behaviour of light in a finite chain of alternat¬ 
ing waveguides with gain and loss have been addressed 
in [18,48]. 


4. Solitons in couplers and periodic 
potentials 

Since the discovery of solitary waves in shallow water by 
Scott Russell in 1834, solitons have attracted a lot of atten¬ 
tion due to their importance in many nonlinear physical sys¬ 
tems, ranging from solids to Bose-Einstein condensates [60]. 
In optics, solitons were investigated and described in many 
research papers and books (see, e.g., Ref. [13]). Conserva¬ 
tive solitons requiring balance of nonlinear response and 
medium dispersion usually form families with different am¬ 
plitudes. Nonlinear dissipative systems, however, demand an 
additional balance between gain and loss to support soliton 


Copyright line will be provided by the publisher 







LASER 

&PHOTONICS 

REVIEWS 


12 


Sergey V. Suchkov et al.: Nonlinear PT-symmetric photonic systems 



Figure 13 (a) Schematic of a planar PT-symmetric coupler and 

(b) a cross-section of its discrete analog. Red color indicates gain 
(black circles in (b)), while blue color indicates loss (white circles, 
respectively). In (b) couplings between waveguides are marked 
as C and k. Adopted from [59,65]. 



Figure 14 Unstable dynamics of a symmetric soliton with (a) 
transformation to a breather, and (b) PT-symmetry breaking (blow 
up). The u component is shown. Adopted from [63]. 


solutions [14,15]. This latter requirement is usually satisfied 
only for particular soliton amplitudes and shapes, thus no 
continuous families can generally be found. In their turn 
PT-symmetric systems belong to a subclass of dissipative 
systems, however due to the symmetry property they can 
commonly support families of solitons with different ampli¬ 
tudes [61]. Accordingly, PT symmetric systems are unique, 
combining certain features of conservative and dissipative 
systems. 

In this section, we discuss how PT-symmetry affects the 
soliton properties compared to conservative or dissipative 
systems. 


4.1. Solitons in PT-symmetric couplers 


First, we consider a planar PT-symmetric coupler where 
light experiences diffraction in the x direction and propa¬ 
gates along the z coordinate [59,62-64]. Figure 13 (a) show 
a schematic pretention of the coupler where the red plane 
stands for a waveguide with gain, while the blue plane shows 
a lossy waveguide. 

This structure can be described by the normalized sys¬ 
tem of coupled-mode equations: 

iu z + u xx + | u | u — iyu + v = 0, (29) 

iv z + v xx + |v| v-\-iyv-\-u = 0, 

where u(z,x) and v(z,x) are the mode amplitudes in the 
waveguides with gain and loss, respectively, y > 0 is the 
gain/loss coefficient. System (29) supports two types of 
soliton solutions: symmetric and antisymmetric [62, 63], 
termed by analogy with conservative systems (when y = 
0) [13]. These two solitons can be presented in the following 
form: 


( 



= Asech 


A 

L7! 


(x — Vz) 





(30) 


Here A, V, j8±, are the soliton amplitude, velocity and prop¬ 
agation constant, respectively, satisfying the relation f3± = 

—A 2 — cos8 ± +V 2 /4+A 2 /2 with cos<5± = =b \f \1 — y 2 and 
sin<5± = — y. The signs =b correspond to the symmetric 
and antisymmetric solitons, respectively. Stability anal¬ 
ysis reveals that symmetric solitons are stable if A 2 < 


(4/3) yj\ — y 2 , while antisymmetric solitons are always un¬ 
stable. However, their lifetimes are exponentially long if 
amplitudes are small enough, so they could be considered to 
be stable in practice. Generally speaking, it is easy to show 
that under assumption u[z) = v(z) exp(/<5±) Eqs. (30) can 
be reduced to a single conservative nonlinear Schrodinger 
equation. It means that, if the solitons (30) are stable, then 
their properties do not differ from ones of their conserva¬ 
tive counterparts, and full energy of the system is exactly 
conserved. However, unstable behaviour of such solutions 
may demonstrate unique features. If the soliton amplitude 
exceeds the critical value, two types of unstable dynam¬ 
ics for each type of solitons are possible: the emergence 
of breathers with oscillating energy or unbounded growth 
of the mode intensity in the waveguide with gain. In the 
conservative case, the latter type of dynamics leads to for¬ 
mation of an asymmetric soliton, which concentrates mostly 
in one waveguide. Since the asymmetric solitons break a 
balance of gain and loss in PT systems, it leads to the PT- 
symmetry breaking. Examples of these unstable behaviors 
are presented in Fig. 14 for a symmetric soliton. 

The considered system supports not only fundamental 
solitons, but also stable A-soliton complexes as well [64]. 
It turns out that unlike the fundamental solitons considered 
above, the stability region for the antisymmetric two-soliton 
is considerably larger than for its symmetric counterpart. 

As mentioned above, unstable solitons can transform 
into breather-like objects [59]. Such nonlinear modes are 
stable on the distances of z < £~ 2 , where e 1 / 2 gives the 
scale of the amplitude of the small-amplitude breathers. 
The breathers are more common objects in the planar PT- 
symmetric coupler in the following sense. Propagation of 
low intensity light (determined by £) is governed by a con¬ 
servative system of amplitude equations obtained in [59]. 
This system admits soliton solutions, which correspond to 
breather solutions in original model (29) as well as degener¬ 
ate soliton solutions, which are symmetric and antisymmet¬ 
ric solitons of the model (29). Breathers can be produced 
resulting from collision of symmetric and antisymmetric 
solitons [see Fig. 15(a)] as well as at fission of unstable soli- 
tons[see Fig. 14 (a)]. Remarkably, although the breather en¬ 
ergy is not conserved during propagation, it is conserved on 
average. An interesting example of breather dynamics with 
merging of two breathers into one is shown in Fig. 15(b). 
Breathers dynamics is very sensitive to a relative phase of 
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Figure 15 (a) Breather emergence in collision of symmetric and 
antisymmetric solitons. (b) Fusion of two breathers into one. Only 
u component is shown. Adopted from [59]. 



Figure 16 Propagation of ’’Zeno” soliton in a coupler composed 
of lossy and conservative waveguides. Input signal was injected 
into the conservative waveguide only (u component), (a) Loss 
coefficient y = 1, in (b) y = 10. Adopted from Ref. [67]. 

interacting breathers [66]. Thus, for fixed initial conditions 
varying only a relative phase of the colliding breathers the 
different types of dynamics such as merging two breather 
into one (temporal or permanent), elastic collision, colli¬ 
sion with energy exchange and PT-symmetry breaking are 
observed. 

A nonlinear dissipative coupler can demonstrate unique 
features due to PT-symmetric behaviour. ”Zeno”-like effect 
was observed in numerical simulation [67] for a system ob¬ 
tained from Eq. (29) by replacing the waveguide with gain 
by conservative one. The Zeno phenomenon, introduced 
in quantum mechanics [68] and consisting in strong sup¬ 
pression of the decay of an unstable particle by means of 
frequent measurements. In terms of considered model more 
frequent measures correspond to stronger loss in the system. 

Generally speaking, the new system is no longer PT- 
symmetric, however, a linear analogue of such a system 
can be brought to PT-symmetric view by simple change of 
variables [7]. As discussed above, if the gain/loss parameter 
is below the PT-symmetry breaking threshold there appears 
beating of light between the waveguides, while in the oppo¬ 
site case this beating is suppressed, and light gets localized 
in the waveguide with gain. At the same time, the systems 
with loss only also exhibit similar dynamics. Thus, stronger 
loss in the system suppresses beating of light, locking it in 
the conservative waveguide and making the whole system 
more transparent (see Fig. 16). 


A discrete analogue of the model (29) has been consid¬ 
ered in Ref. [65]. The structure is a chain of linearly coupled 
PT-symmetric dimers shown schematically in Fig. 13(b). An 
important distinction of the discrete model from its continu¬ 
ous counterpart is that the former model possesses discrete 
linear spectrum providing additional possibility for solitons 
to be stable. Two types of solitons, symmetric and antisym¬ 
metric, were constructed numerically, and their properties 
were investigated. It turns out that the discrete PT-symmetric 
lattice supports freely moving solitons. Solitons behave 
quite similar to solitons in conservative systems provided 
they are stable and their amplitude does not exceed a critical 

value determined as |A cr i t | 2 = \J 1 — y 2 . If soliton amplitude 
exceeds A cr i t , the u component of the soliton, which experi¬ 
ences gain, grows exponentially, while the component with 
loss, v, quickly decays. 

The discussed model (29) as well its discrete counterpart 
do not admit dark soliton solutions [see for example [65]]. 
However nonlinear coupling between waveguides (cross¬ 
phase modulation) can facilitate the formation of dark and 
bright solitons. Also stable rogue waves of the Peregrine 
type were identified [69,70]. It was shown that gain/loss pa¬ 
rameter modifies the region of stability of both background 
of dark solitons as well as solitons themselves. 


4.2. Solitons in localized potentials 

As was mentioned above, the necessary condition for one¬ 
dimensional optical system to be PT symmetric is U (x) = 
£/*(—. x) [71]. The presence of nonlinear response of a 
medium, generally speaking, changes the real part of system 
potential making eigenvalues of the Hamiltonian complex¬ 
valued. However, judicious distribution of an intensity in 
such systems may retain eigenspectrum to be real. 

Fight propagation through the waveguide with a mod¬ 
ulated refractive index in transverse direction (x) and ex¬ 
periencing gain and loss is described by the normalized 
equation 

iVz + Vxx + U(x)\i/+o\y\ 2 y = 0, (31) 

where U(x) = V (x) + iW (x) satisfies the condition above 
for PT-symmetric systems. In this section we consider the 
case, when both V (x) and W (x) are localized in transverse 
direction v. 

Stable solitons in PT-symmetric systems with distributed 
gain and loss were first reported in Ref. [71]. Particu¬ 
larly, the potential of Scarff form V = Vo sech 2 (v), W (x) = 
Wo sech(v) tanh(v) has been considered for focusing type of 
nonlinearity [see Fig. 17 (a,b)]. It has been revealed that 
the system possesses real linear spectrum if Wo < Vo + 0.25. 
However, even if this condition is broken (linear spectrum 
is complex), nonlinear modes can alter the amplitude of re¬ 
fractive index due to nonlinearity, thus bringing the system 
into PT-symmetric phase. Some analytical solutions have 
been obtained in [72] for ID and 2D PT-symmetric wells 
and lattices. Although the soliton in Fig. 17(b) looks similar 
to ones determined by Eq. (30) and depicted in Fig. 14, in 
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Figure 17 Solitons in localized PT-symmetric potentials, (a) The 
Scarff type potential: blue indicates a real part of the potential, red 
- its imaginary part, (b) Soliton profile and its intensity distribution 
(adopted from [71]). (c) Parabolic trap potential considered in 
Ref. [73] with £2 = 0.1 (blue) for two cases of gain/loss strength: 

£ = = 0.064 (red) and £ = 0.64 (green), (d) Dark-soliton 

dynamics corresponding to £ = 0.064 (top) and £ = 0.64 (bottom). 
Adopted from [73]. 




Figure 18 (a) Soliton profile (blue curve - real part, red curve 

- imaginary part) and potential (refractive index) distribution, (b) 
Stable soliton propagation in a PT-symmetric periodic potential. 
Adopted from [71]. 

sidered case, then instead of the exponential growth one 
observes the soliton generation from a background state 

[see Fig. 17(d), bottom panel]. This critical point £ crit cor¬ 
responds to blue-sky bifurcation of the background state 
and the fundamental soliton branch, after which they cease 
to exist. The same type of bifurcation is also observed for 
soliton branches of higher orders, which are found for larger 
values of £, i.e. 2-soliton branch bifurcates with 3-soliton 

branch at some £ crit > £ crit etc. Thus, for larger values of £ 
(strong gain/loss) we may observe spontaneous generation 
of multiple dark solitons from an input beam. 


4.3. Solitons in periodic potentials 


fact they are completely different. Each of the components, 
u and v, of the soliton (30) experiences pure gain or loss 
respectively, and the energy flow is directed from the gain 
waveguide to the lossy one. The soliton presented in Fig. 17 
has only one component, which is partially situated in both 
gain and loss regions. This soliton has the transverse energy 
flow, which strongly affects its properties, e.g. stability and 
mobility. 

In the case of defocusing nonlinearity [74-76], gray 
solitons, nonlinear guided modes as well as bright solitons 
in nonlocal nonlinear media, have been found to be stable 
in PT-symmetric localized potential of the Scarff type. 

In Ref. [77] stable fundamental soliton, as well as 
two- and three-soliton solutions, have been reported for 
the Gaussian PT-symmetric potential of the form U(x) = 

(1 + iWox ) exp (—v 2 ). 

An interesting study of the behaviour of dark solitons 
near the phase transition points was reported in Ref. [73]. 
A parabolic trap potential of the form V(x) = l/2£2 2 v 2 
and W(x) = £vexp(— x 2 /2) was considered [see Fig. 17(c)] 
and it has been found that dark solitons are stable if 
the amplitude of the imaginary potential, £, is less than 

e crit = a/ 5/I2Q. An unstable soliton of an oscillating 

type is shown in Fig. 17(d) (top panel) for £ = £ c /| and 
O, = 0.1. Above this critical value, the dark-soliton back¬ 
ground demonstrates exponential growth corresponding to 

the blow-up due to PT-symmetry breaking. However, if £ 

( 2 ) 

exceeds some critical value, namely £ crit = 0.62 in the con¬ 


Next, we consider soliton solutions in structures having peri¬ 
odic modulation of both a real part of the refractive index as 
well as gain/loss regions. The governing equation for light 
propagation has the form of Eq. (31), but potential U (x) is 
periodic. Linear properties of such a system have been inves¬ 
tigated in Ref. [78]. It has been shown that nonorthogonality 
of thed Floquet-Bloch modes for PT-symmetric potentials 
leads to effects such as double refraction, power oscilla¬ 
tion (due to unfolding of nonorthogonalized Floquet-Bloch 
modes), nonreciprocity, and secondary emission in 2D case. 
These linear attributes of PT-symmetric system may also 
be observed in nonlinear models and are fundamental phe¬ 
nomenons that distinguish PT-symmetric systems from con¬ 
servative ones. 

In the case of self-focusing nonlinearity, <7 = 1, soli¬ 
ton solutions have been investigated for the potential of 
the form V(x) = cos 2 (v), W(x) = Wbsin(2x) [71]. It is in¬ 
teresting that nonlinear response in such systems can shift 
PT-symmetry breaking threshold enabling one to find non¬ 
linear mode with real propagation constant even though the 
associated linear spectrum is complex. It turns out that the 
solitons can be stable over a wide range of model param¬ 
eters. Furthermore, narrower self-trapped waves are more 
stable and it is possible to find stationary modes even above 
PT-breaking threshold, but this class of solitons unfortu¬ 
nately is unstable. Potential distribution as well as stable 
soliton propagation are shown in Fig. 18. 

A system of superlattices with V (x) = £ sin 2 [x + n/2] + 
(1 — £)sin 2 [2(v + 7t/2)\ and W(x) = Wbsin(2v) has been 
studied in [79]. Parameter £ here controls the relative 
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Figure 19 Nonlinearity induced transition from PT broken to full 
PT phase, (a-c) Invariant propagation of signal, the initial intensity 
of which is high enough to transit the system to PT symmetry 
phase, (a) Packet intensity, (b) power, (c) center of mass, (d-f) The 
same as in (a-c) but for the packet with insufficient initial intensity 
for transition. Adopted from [81]. 


strength of the superlattices. Stable solitons propagating 
in the low power region were found in semi-infinite gap. 
The relative strength of the superlattices may significantly 
vary PT breaking threshold as well as region of soliton 
stability. 

For the defocusing nonlinearity, a = — 1, multi-peak 
gap solitons exist and are stable over a wide range of 
parameters in a first gap of the potential of the form 
Vb sin 2 (x) H-Wo sin(2jc) [80]. It is remarkable that stable 
solitons with larger number of peaks exist for a wider range 
of propagation constant. 

Interesting finding has been obtained in Ref. [81]. It 
turns out that nonlinearity may cause a transition from bro¬ 
ken to PT symmetric phase for focusing nonlinearity and 
vice versa for defocusing one. Nonlinear band structure has 
been obtained for different strength of nonlinearity and it 
has been revealed, that the greater the nonlinearity strength, 
the smaller region in the Brilliouin zone where the eigenen- 
ergies are complex. When power reaches a certain threshold, 
the whole spectrum becomes real. Thus, two different dy¬ 
namical regimes for nonlinear wave packet near the phase 
transition point were observed. First, if an initial power of 
a propagating packet is large enough for a transition of the 
system to PT-symmetric phase, then intensity and power are 
invariant throughout propagation. In Figs. 19 (a-c) an inten¬ 
sity, power and center of mass of the propagating packet are 
presented for this case. Second, if the initial power of the 
packet is too low to induce transition, then the intensity and 
power start to grow [see Fig. 19 (d-f)]. After power passes 
the level at which spectrum becomes real it continues to 
grow until it reaches some maximum value (indicated by 
black line in Figs. 19 (e,f)) and after that it starts to decay 


bringing the system to the original state. This dynamics can 
be explained as follows: the evolution starts an the center 
of unit cell (we consider here only one period of the poten¬ 
tial. Since the initial power is low (linear regime) we can 
decompose the packet into two eigenmodes, one of which is 
growing and the second one is decaying. The growing eigen¬ 
mode resides mostly in the gain half of the unit cell, while 
decaying - in the loss half. After the decaying mode dies out, 
the center of mass of the packet is shifted along the potential 
minima to the gain region. Then it starts to move in opposite 
direction along with the power growth and ends up in the 
loss region. The power of the packet starts to decay, bringing 
the system to PT broken phase and the process repeats again. 
Both of these dynamical regimes are stable and numerical 
simulations indicate that such propagation should continue 
indefinitely. The second type of dynamics is associated with 
any continuous physical systems and is not observed in sys¬ 
tems modeled by a standard tight-binding model, since for 
nonlinearity induced transition the space overlap of modes 
corresponding to gain and loss should be presented. This 
effect allows to design systems that can switch from being 
power conservative to amplification regime. 

Behaviour of wavepackets near the phase transition 
point has been studied in [82]. In this case the envelope 
dynamics of the propagating packets can be described by 
nonlinear Klein-Gordon equation (Eq. (13) in [82]). Based 
on this envelope equation a variety of phenomena such as 
wave blow up, periodic bound states and solitary waves are 
predicted and corroborated by comparison with numeric 
simulation of full model (31). 

A discrete counterpart of the PT-symmetric periodic 
potential can be realized by appropriate alteration of cou¬ 
pled waveguides with gain and loss. Such discrete structures 
are easier to fabrication then the waveguides with continu¬ 
ously varying refractive index as well as gain/loss strength. 
Solitons in discrete PT-symmetric lattices have been investi¬ 
gated in a series of works [18, 83, 84]. A chain of alternating 
PT-symmetric couplers was considered in [83]. The field 
evolution along the waveguides is described by the system 


duj 

dz 

dvj 

l ~dz 


T ijUj -f-Cl vj -f- C2V j— 1 -f- 1 Uj\^Uj — 0, 
- iYVj+C\Uj+C2Vj+i + I Vy | 2 vy = 0. 


(32) 


Here uj and vj are the mode amplitudes in waveguides with 
loss and gain, respectively, C\ and C2 - coupling coefficients. 
In the linear regime the system has real eigenspectrum if 
||Ci|-|C 2 || > y|. It has been shown numerically that this 
lattice supports stable soliton solutions. These solitons might 
be actively reshaped or switched to PT-broken phase by a 
small-amplitude modulation of the gain/loss parameter y 
along the propagation direction z. Later, in the work [84] 
an analytical proof of existence of such solitons was given 
and classification of solution families and their stability was 
provided. Moreover, by analytical continuation from the 
anticontinuum limit it is possible to show that localized 
discrete solitons exist in a network of PT-symmetric clusters 
(oligomers). 
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4.4. Generalized nonlinear potentials 

In this section we consider PT-symmetric systems with 
generalized nonlinear response such as nonlocal nonlinear¬ 
ity, cubic-quintic competing nonlinearity, and mixed linear- 
nonlinear PT-symmetric potentials etc. The held evolution in 
one-dimensional systems with a general type of nonlinearity 
can be described by the following normalized equation 

i-^- + -+ [V (x) + iW (x)] y/ + G(|i/a| 2 ,x)i/a = 0, (33) 

where G(|i/a| 2 ,x) is a functional determining the nature of 
nonlinear response. 

Nonlocal nonlinear response appears in such systems 
as photorefractive crystals, nematic liquid crystals, lead 
glasses, etc. This nonlinear response can significantly mod¬ 
ify the soliton properties such as existence, stability, and 
mobility. The function G satisfies the following equation 
G-d(d 2 G)/(dx 2 ) = Iv'l 2 . Here d is a degree of nonlocal¬ 
ity of the nonlinear response. Several papers were devoted to 
the study of gap solitons and defect modes in self-focusing 
nonlocal nonlinear media. In Ref. [85], the existence and 
stability of defect modes were studied for different values 
of the nonlocality parameter d , as well as for different types 
of defects. In Ref. [86], gap solitons were studied in dual 
periodic lattices, while Ref. [87] was devoted to the study 
of defect modes in PT-symmetric superlattices. In Ref. [88] 
the effect of spatially modulated nonlocal nonlinearity has 
been addressed. It was shown that there exist regions where 
gap solitons and defect modes are stable, and that the degree 
of nonlocality affects this region drastically. 

Defocusing type of nonlinearity has been considered in 
[89]. The existence and stability of gap solitons in media 
with diffusive-like nonlinear response were determined. The 
gap solitons are always oscillatory unstable for any degree 
of nonlocality in the presence of an imaginary potential. 

It has also been shown that PT-symmetric periodic lat¬ 
tices can support a variety of stable vector solitons consist¬ 
ing of multiple bright beams in both focusing and defocus¬ 
ing media [90] . Dynamics of such solitons is described by 
Eq. (31) where \jf = (VO; V^) is a two component vector. 
In the component form Eq. (31) becomes a system of two 
incoherently coupled-mode equations. Stable solitons can 
be found even if both components have different propaga¬ 
tion constants, but to satisfy the PT-symmetry condition 
these solitons must have equal number of humps in each 
component. It also may happen, that one of the compo¬ 
nent is unstable itself, while the other is stable. In this case 
the soliton can be stable due to the coupling of these two 
components. Analysis of an effect of gain/loss strength on 
soliton stability revealed that stronger gain/loss coefficient 
can destabilize soliton in focusing media and stabilize it in 
defocusing one. 

The case of spatial periodic modulation with Kerr- 
type nonlinearity also including gain/loss modulation has 
been considered in Ref. [91]. In this case G(|i/a| 2 ,x) = 
V[(x) + iW (x)] | xj/\ 2 with zero linear potential and the func¬ 
tions V (x) and W (x) have the same modulation period. Anal¬ 
ysis of the soliton existence and stability reveals that in the 





Z 


Figure 20 Bragg solitons in PT-symmetric potential. (a,b) For¬ 
ward and backward waves below PT-symmetry breaking threshold, 
(c,d) the same but at the PT-symmetry transition point. Initially 
only forward wave was excited. Adopted from [96]. 


case of zero real part of the nonlinear potential there exists 
a non-vanishing region of parameters where stable solitons 
exist. The modes with the width less then a half of the lat¬ 
tice spacing appear to be stable. The system supports stable 
multihump solitons. Stable solitons were found in the lattice 
with mixed linear-nonlinear PT-symmetric optical poten¬ 
tials [92]. Some exact solutions to the system with mixed 
linear-nonlinear PT-symmetric potentials were obtained in 
Ref. [93]. 

An effect of multistability has been observed in a 
model with competing cubic-quintic nonlinear response 
with G = (|i/a| 2 — |i/a| 4 ) [94,95]. At this type of nonlinear 
response for low intensities of light the focusing nonlinear¬ 
ity prevails on defocusing one, while for high intensities 
defocusing dominates. It has been shown [95] that solitons 
with different symmetries and different number of humps 
can propagate with the same propagation constant. Another 
interesting point is that both even and odd hump solitons 
can be stable over some range of propagational constants, 
which is in contrast to conservative case, where only lowest 
branches of solitons of both types can be stable. Thus imag¬ 
inary part of the potential can stabilize solitons in the first 
gap. 

In the work [96] the effect of PT-symmetric potential 
on Bragg solitons has been studied. These solitons are com¬ 
posed of interlocked both forward and backward propaga¬ 
tion modes. Such interlocking open up forbidden band gap 
allowing energy transport. It is interesting that at the exact 
PT phase transition point the exchange of energy between 
forward and backward modes ceases, and full energy is con¬ 
served. Below PT breaking threshold, there is oscillatory 
power exchange between forward and backward modes [see 
Fig. 20]. 

Interesting observation has been made in [61]. It has 
been shown analytically and corroborated numerically that 
for ID nonlinear Schrodinger equation PT-symmetry is a 
necessary condition for existence of a continuous family of 
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solitons bifurcating out from real linear spectrum. In con¬ 
trast for other dissipative systems based on non-PT symmet¬ 
ric potentials with all-real spectrum only isolated solitons 
can exist. From the other hand, the nonlinear Schrodinger 
equation with a complex PT potential cannot admit a con¬ 
tinuous family of non-PT-symmetric solitary waves [97]. 

4.5. Two-dimensional potentials 

Here we describe the effect of PT-symmetric potentials on 
nonlinear modes in 2D case. Due to an additional degree of 
freedom there appear new types of nonlinear modes such as 
spatiotemporal solitons and vortex solitons. 

Propagation of light in 2D nonlinear PT-symmetric sys¬ 
tems with Kerr-type nonlinear response can be described by 
the equation similar to Eq. (31) 

i\l/ z + V 2 xi/+U(x,y)\i/+a\\i/\ 2 xi/ = 0, (34) 

with two dimensional PT-symmetric potential U(x,y) = 
V(x,y)+iW(x,y). 

One of the first examples of 2D solitons in a PT- 
symmetric periodic potential has been provided in [71]. 
Later, in Ref. [98] the multipeak gap solitons in a PT- 
symmetric periodic potential were considered. It was shown 
that the system (34) with defocusing Kerr-type nonlin¬ 
earity and PT-symmetric potential of the form V(x,y) = 
Vb{cos 2 (x) + cos 2 (y) + /Wo[sin(2x) + sin(2y)]} admits sta¬ 
ble multipeak gap solitons. Remarkably, that not only even 
peak soliton might be stable, but odd-peak solitons as well. 

The system (34) possesses distinct dynamics near the 
PT phase transition point. For the potential U(x,y) = 
V 0 2 [cos(2v) + cos(2y) + iWo( sin(2v) + sin(2y))] it has been 
shown, that wave packets experience pyramid like diffrac¬ 
tion near the transition point Wo = 1 [99]. Near this point 
diffraction surface intersects itself four-fold and wavepacket 
dynamics can be described by fourth-order envelope equa¬ 
tion [see Eq. (11) in Ref. [99]]. It was found that below 
certain threshold (which is equal to 3.2 in the considered 
case) a wave packet exhibits discrete (pyramid type) diffrac¬ 
tion as in the linear case. Above this threshold an exponen¬ 
tial growth is observed. It is quite surprising that the wave 
packet dynamics is insensitive to the type of nonlinearity 
(compared to one-dimensional case considered in Ref. [82]) 
as follows from the middle and right columns of Fig. 21. In 
Fig. 21, nonlinear dynamics for different initial amplitudes 
A and two types of nonlinearity is presented. Upper row is 
obtained from simulations of the envelope equation, lower 
row from full model simulations. 

The PT-symmetry potentials can be taken as a base for 
systems, which, generally speaking, are not PT invariant. In 
the work [100] the Anderson localization of light has been 
studied in the system (34) with PT symmetric potential, 
and for introduced disordered potential depth Vo along the 
waveguide. Although there are exponentially growing (in 
z direction) modes in the system, for low level of disorder 
as well as for short propagation distance it is possible to 
observe an effect of PT-potential on transverse localization 
of light. It turns out that PT-symmetric potential enhances 


=2 Aq =3.3 Aq —4 


Initial Wavcpackct Focusing (a = 1) Defocusing (a = -1) 


Figure 21 Pyramid like diffraction in 2D PT-symmetric potential 
near the transition point Wq = 1. Left column - below transition 
point, middle and right above. Upper row corresponds simulation 
of the envelop equation, lower row - full model. Adopted from [99]. 

the transverse localization of light in comparison with con¬ 
servative systems as well as systems with just gain or loss. 

One promising observation that could find some ap¬ 
plications was suggested in Ref. [101]. A partially PT- 
symmetric potential, satisfying the following property 
V*(x,y) = V(— x,y), can possesses all-real spectra and sup¬ 
port continuous families of stable solitons as in case of full 
PT-symmetry. Remarkably, this soliton family may exhibit 
multiple power branches with upper branches more stable 
then the lower ones. These results expand the concept of PT 
symmetry in multidimensions. 

By analogy to the ID discrete potential considered 
in [65], a 2D discrete nonlinear PT-symmetric lattice has 
been considered in [102]. An array of coupled PT-symmetric 
dimers can be realized in form of a set of dual-core wave¬ 
guides embedded into a photonic crystal. This lattice sup¬ 
ports stable symmetric and antisymmetric fundamental and 
vortical discrete solitons. A symmetric fundamental soliton 
represents the ground state of the system, and this soliton 
can be stable at lower values of the total power, while anti¬ 
symmetric fundamental solitons as well as on-site-centered 
vortices tend to be stable at higher powers. 

5. Scattering on PT-symmetric potentials 

Various optical systems can be probed in terms of light 
scattering, which delivers important information about their 
spectral properties and reveals the consequences of their 
symmetries. PT potentials enable new regimes of linear scat¬ 
tering, which cannot be achieved in conservative structures, 
and furthermore facilitate efficient light switching including 
nonreciprocal response in the nonlinear regime. 

5.1. Linear scattering and invisibility 

Combination of active and lossy optical elements can be 
designed in a special way to achieve counterintuitive effects 
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including absence of reflections and invisibility [42, 103, 
104]. The scattering properties of a one-dimensional system 
can be characterized by the matrix 


S = 


T R r 
R l T 


(35) 


where T and R r l respectively are the (complex) transmission 
and right/left reflection amplitudes, so that the transmission 
and right/left reflection coefficients are given by \T \ 2 and 
\R r,l \ 2 . The scattering matrix is formulated taking into ac¬ 
count equality T l = T r = T, which reflects the reciprocity 
of transmission for any linear, stationary, and non-magnetic 
medium. 

The scatterer is called reflectionless from the left (right), 
if R l = 0 and R 1 fz 0 (R 1 = 0 and R l 0). The scatterer is 
called invisible from the left (right), if it is reflectionless from 
the left (right) and in addition T — 1. A thorough assessment 
of the role of PT-symmetry in the phenomenon of invisibil¬ 
ity of scattering potentials in one dimension has been done 
by Mostafazadeh [103]. PT-symmetric nature of invisibility 
is established in this work by the following invisibility theo¬ 
rem. Consider a general real or complex one-dimensional 
scattering potential v that is defined on R. Let v PT be the 
PT transform of v that is given by v PT (z) := v(— z)*, and k* 
be a positive real number. Then for an incident wave with 
wavenumber k, the following equivalent statements hold: 

(i) v is invisible from the left (or right) for k = k* if and 
only if so is v PT . 

(ii) v is invisible from the left (right) for k = k* if and 
only if v* is invisible from the right (left) for k = k*. 

Invisibility does not require PT-symmetry of v and for 
non-PT-symmetric potentials this theorem implies a pairing 
of the left/right-invisible configurations that are related by 
the PT transformation. Interestingly, the statement of the 
invisibility theorem also holds, if the term ’’invisible” is 
replaced with ’’reflectionless” [103]. 

A simple two-layer structure shown in Fig. 22 represents 
a basic, exactly solvable, and experimentally realizable sys¬ 
tem which can realize invisibility [103]. Infinite planar slab 
of optically active material consists of two layers having 
complex refractive indices n\ and n 2 and equal thickness 
of L/2. The wave equation for a linearly polarized time- 
harmonic electromagnetic wave that propagates along the 
direction normal to the slab reads 

'i>"(z)+k 2 n{z) 2 'i>(z) = 0, (36) 

where 

( n\ for — L/2 < z < 0, 
n(z) = < n 2 for 0<z<L/2, (37) 

l 1 for \z\ > L/2. 

The amplitudes of plane waves depicted in Fig. 22 (a) are 
related through a 2 x 2 matrix M 



Figure 22 Linear scattering regimes, (a) Schematic of scat¬ 
tering by active two-layer slab, as formulated in Eqs. (37), (38). 
(b-d) Experimental demonstration of unidirectional reflectionless 
scattering by active Bragg gratings, (b) Design of passive structure 
and SEM picture of the fabricated device, (c) Measured spectrum 
of contrast ratio of reflectivities for left and right incidence, (d) Sim¬ 
ulated electric field distribution for incidence from the left (top) and 
right (bottom). Adopted from Ref. [104]. 


Particular cases of left-incident and right-incident scattering 
solutions can be written as [105] 



and 


ci(e lkz + rie lkz ) for z<—L/ 2, 
citje lkz for z>L/ 2, 


(39) 


Viz) 


c t p i k ^ 

L plpC' 

c r (e~ lkz -f- r r e lkz ) 


for z < -L/2 , 
for z > L/2 , 


(40) 


where the incident, transmitted and reflected light intensities 
are = \c ^ r | 2 , T Lr = \t^ r \ 2 , and R Lr = \r^ r \ 2 , respectively. 
The (complex) transmission and left/right reflection ampli¬ 
tudes can be determined by the elements of matrix M (note 
that det(M) = 1) as follows 


1 M 21 M\2 

ti=t r = -, n =-, r r =- 

M22 M22 M22 


(41) 


with Mij defined by Eq. (16) in Ref. [103]. A spectral sin¬ 
gularity appears whenever one can satisfy M 22 = 0 for a 
real k, then the transmission and reflection amplitudes di¬ 
verge. Then, the necessary and sufficient conditions for the 
unidirectional invisibility are found as 


n+ cos a + — n 2 _ cos a _ = ( n\ — n 2 _) cos (Lk ), 

r\ r\ 

h+n+ sin +«_«_ sin a- = (ni J_ — n_) sin(Lfc), 

cosa + ^cosa_, (42) 


where 


n±=n\±n 2 , a±=n±kL/ 2, h±=n\n 2 ±l- (43) 


= M 


(38) 


The two-layer slab is invisible from both directions provided 
that n\ and ^2 are rational numbers, and the wavelength A 
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at which the device is invisible is such that L is the half 
integer multiple of A [103]. Moreover, existence of invisible 
configurations that remain practically reflectionless within 
a very wide spectral range was demonstrated [103]. 

Unidirectional invisibility can be also achieved for PT- 
symmetric locally periodic potentials describing Bragg grat¬ 
ings with active layers [42] 


v(z) 


(do + ae 2, P z for \z\<L/2, 
\0 for \z\>L/2, 


(44) 


where (Xo, oc, jS > 0, and L > 0 are particular real parameters. 
It was then shown that invisibility can be achieved only 
for structure length below a certain threshold [106], and 
analytical solutions were derived in Ref. [107]. 

The invisibility phenomenon was demonstrated exper¬ 
imentally in Ref. [104]. The optical non-Hermitian parity¬ 
time system was fabricated with only absorptive media on 
the Si-on-insulator (SOI) platform, see Fig. 22(b). This struc¬ 
ture effectively realized a parity-time-symmetric periodic 
potential according to Eq. (44), but in a lossy background. 
Accordingly, the transmission is attenuated, but unidirec¬ 
tional reflection can still be observed. Indeed, experimen¬ 
tal measurements show large contrast for reflections for 
incidence from the left and right, see Fig. 22(c). The dif¬ 
ference occurs due to distinct distributions of the electric 
held, which are visualized based on numerical modelling in 
Fig. 22(d). 

Suchkov et al. have studied wave scattering on a do¬ 
main wall introduced into a linear array composed of PT- 
symmetric dimers by switching the gain and loss in a half 
of the array [108]. The linear spectrum of the defect-free 
array includes high-frequency and low-frequency branches. 
The two major effects have been revealed: amplification of 
both reflected and transmitted waves and excitation of the 
reflected and transmitted low-frequency and high-frequency 
waves by the incident high-frequency and low-frequency 
waves, respectively. 

General approach to scattering theory for optical PT- 
symmetric systems with a wide range of geometries beyond 
the one-dimensional setting has been discussed in a recent 
work by Schomerus [109]. In particular, it is demonstrated 
that in the presence of the magneto-optical effects the sym¬ 
metry of optical device becomes important. In addition to 
the PT-symmetry, the PTT'-symmetry is introduced to ac¬ 
count for the variants of the geometrical and time-reversal 
operations. The latter symmetry imposes the same spectral 
constraints as the former one but, generally speaking, a dif¬ 
ferent condition on the magnetic held. Scattering theory is 
applied to study spectral features of the open optical sys¬ 
tems possessing different symmetries. This leads to effective 
models that can be formulated in the energy domain (via 
Hamiltonians) and in the time domain (via time evolution 
operators) giving the information about the related energy 
and time scales. 


5.2. Nonlinear nonreciprocity 

Whereas transmission is known to be reciprocal in linear 
systems as discussed in the previous section, nonlinearity 
can lead to nonreciprocal transmission which depends on the 
direction of incidence, i.e. 7) y T r . Nonlinear systems can 
realize optical isolation by predominantly allowing only one¬ 
way transmission, i.e. it allows light to pass through, say, for 
left incidence, whereas it is blocked for light coming from 
the other side. Such nonlinear optical diode have been exten¬ 
sively investigated in conservative optical systems, however 
active PT structures can be used to dramatically reduce the 
optical power threshold for nonlinear nonreciprocity [110]. 

The effect of nonreciprocity was first observed experi¬ 
mentally for electric current transmission lines attached to 
PT-symmetric electronic circuit resonators playing the role 
of nonlinear scatterers [111, 112]. As a reference model, 
Bender et al. used coupled nonlinear electronic Van der Pol 
oscillators with anharmonic parts consisting of a comple¬ 
mentary amplifier (gain) and a dissipative conductor (loss) 
combined to preserve PT-symmetry. A striking feature ob¬ 
served in measurements was that the transmittance from 
left to right differs from the transmittance from right to left. 
The phenomenon is most pronounced in the regions of the 
resonances distorted by the nonlinearity, around 39.5 kHz. 
Furthermore, it was found that increasing the gain or loss 
parameter, maintains or even enhances the transmitted in¬ 
tensities while it leaves unaffected the resonance position. 
This has to be contrasted with other proposals of asym¬ 
metric transport which are based on conservative nonlinear 
schemes (see for example Ref. [113]), where increase of 
asymmetry leads to reduced transmittances. 

Nonreciprocal transmission in optics was observed in 
coupled-resonator structures, where nonlinear effects are 
enhanced due to strong held confinement [110]. The system 
was composed of two coupled microtoroidal whispering- 
gallery-mode resonators [Fig. 23(a)], where one was passive, 
and gain in the other one was provided in the l,550nm wave¬ 
length band by optically pumping Er 3+ ions with a pump 
laser in the 1,460nm band. The coupled-resonator system 
is PT symmetric because under parity reflection P the res¬ 
onators become interchanged and under time reversal T 
loss becomes gain and gain becomes loss. The microtoroids 
were fabricated at the edges of two separate chips placed on 
nanopositioning systems to control precisely the distance 
and accordingly the coupling, see Fig. 23(b). For strong cou¬ 
pling, the modes are approximately PT-symmetric, whereas 
breaking of the mode symmetry is observed when the cou¬ 
pling is decreased beyond a critical value. 

Transmission between ports 1 and 4 was monitored at 
different input power levels. At low powers the input-output 
relation was linear and the system was reciprocal in both the 
broken- and unbroken-symmetry phases. Indeed, a linear 
static dielectric system, even with gain and loss, can not have 
non-reciprocal response. At higher input powers, the system 
remained in the linear regime for the unbroken-symmetry 
phase. However the nonlinear threshold was lower in the 
broken-symmetry phase, due to the stronger held localiza¬ 
tion in the resonator with gain. As a result, nonlinear non- 
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Figure 23 Parity-time-symmetric whispering-gallery microcav¬ 
ity resonators, (a) Schematic of two directly coupled resonators 
and fibre-taper waveguides. Gain is introduced in resonator / 1 R 1 
through Er 3+ doping, (b) Side view of resonators. (c,d) Unidirec¬ 
tional transmission in the nonlinear regime: (c) from port 1 to 4 
and (d) from port 4 to 1. Figures after Ref. [1 1 0]. 

reciprocal effect was observed: the input signal at port 4 was 
transmitted to port 1 at resonance [Fig. 23(d)], but the input 
signal at port 1 could not be transmitted to port 4 [Fig. 23(c)]. 
The advantages of the PT microresonator systems for the re¬ 
alisation of nonreciprocal transmission include a significant 
reduction in the input power to ~ 1 /i W and higher contrast 
with a complete absence of the signal in one direction but 
resonantly enhanced transmission in the other direction. 

PT-symmetric structures based on microresonators can 
also be used to realize nonlinear Fano resonances that may 
give rise to ultralow-power and high-contrast switching 
and non-reciprocity due to their sharp asymmetric line 
shapes [27,114]. Miroshnichenko et al. have considered 
an optical PT-symmetric gain/loss nonlinear dimer coupled 
to a passive chain transmitting linear discrete waves (j> n (t), 
as shown in Fig. 24(a) [27]. Propagation of light in this 
system can be described by the following set of discrete 
nonlinear Schrodinger equations 

i\j/ A =E\i/ A + i( 70 - y> | Wa I 2 ) Ya + Vy/o, (45) 

i% = C(y„- 1 + Wn+\) + VS„ fi (VA + Yb), (46) 

Wb = EYB-i(7o-Yi\YB\ 2 )YB + VYo, (47) 

where the overdot stands for the time derivative, 7o > 0 
accounts for the linear gain and loss acting on complex 
variables 1//4 and i//#, respectively, E is a frequency shift, 72 
accounts for the PT-symmetric nonlinear loss and gain (as 
shown in [27], stable eigenstates are obtained with 72 > 0, 
i.e., if the nonlinear loss competes with the linear gain and 
vice versa), C is the coupling constant in the linear chain, V 
is a coupling coefficient between the dimer and the chain, 
and <5^0 = 1(0) for n = 0(y^ 0). Due to the symmetry, wave 


Figure 24 (a,b) Fano resonances in a linear chain with the 

side-coupled elements featuring the nonlinear PT-symmetry de¬ 
scribed by Eqs. (45-47), adapted from Ref. [27]. (a) Structure 
geometry, the arrows indicate incident, reflected, and transmitted 
waves, (b) Normalized transmission coefficient as the functions of 
input power (Pj = |/| 2 ) for 70 = 0.01, 72 = 0.0001, V = 0 . 2 , C= 1, 
and co = E = 0.1. The nonlinear Fano resonances correspond 
to T = 0. (c) Transmittance and Fano resonances for nonlinear 
microdisks with gain and loss coupled to an access waveguide, 
adapted from Ref. [114]. Dashed lines correspond to transmission 
from the gain side, while solid lines to transmission from a lossy 
side, shown at different input power levels. The transport asym¬ 
metry is as high as 46.5 dB without compromising the outgoing 
optical intensity which is as high as -5 dBs. 

scattering in this system does not depend on the side of 
incidence. 

The solution corresponding to the scattering of incident 
waves with amplitude I on the PT complex is looked for as 

_ ( Ie i(kn-cot) +Re -i(kn+cot) for „ < (y 

<// " “ \ 7yV«-®0 for n> 0 , ’ 

where wave number k > 0 is determined by the dispersion 
equation for the linear chain, k = cos -1 (ft>/2C), while R 
and T are the amplitudes of the reflected and transmitted 
waves. 

In the nonlinear regime with 72 > 0, one can find sym¬ 
metric solutions with 1//4 = —1//5 = —i(/) 9 where (f) is real. 
For (j) 7^ yj 70/72, the symmetric mode exists at co = E, with 
(j) determined by the equation 

y 2 0 3 -y o 0 + V/ = O, (49) 

which yields a single real root for Pj = |/| 2 > 
(4/27)7 q/(V 2 72), and three real solutions (tristability) in 
the opposite case. All these solutions realize the perfect 
transmissivity with T = 1 [the horizontal blue line in Fig. 
24(b)]. 
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In contrast to its linear counterpart, the nonlinear system 
may support complete suppression of the transmission (T = 
0), that is, nonlinear Fano resonances [115]. The nonlinear 
PT-symmetric scatterer attached to the linear chain as shown 
in Fig. 24(a) can support a continuous family of symmetric 
nonlinear Fano resonances with co = E (for the details see 
[27]). The family of the nonlinear Fano resonances with the 
symmetrically excited side-coupled PT elements exists in 
an interval of the intensity of the incident wave Pj, even 
though Fano resonance is usually obtained as an isolated 
solution. 

It was further suggested that nonlinear nonreciprocity 
and optical isolation through Fano resonances can be real¬ 
ized with two nonlinear micro-resonators with gain and loss 
are coupled to a single waveguide [114]. In the nonlinear 
regime, the line-shape and Fano resonance position depends 
on the direction of the incident light, which can provide non¬ 
reciprocal transmission with 47 dBs contrast in the optical 
C-window around 1.55/im optical wavelength. 

Scattering of linear waves propagating in a waveguide 
array with an embedded nonlinear PT-symmetric dimer was 
found to demonstrated distinct features [116]. It was found 
that, under certain conditions, nonlinearity can effectively 
suppress the nonconservative effects, which are generated 
by the linear gain and loss terms. In particular, at high inten¬ 
sities, the incident wave can almost entirely transmit across 
the defect. 

Linear PT-symmetric systems can demonstrate spectral 
singularities leading to infinite transmission and reflection 
coefficients [105]. Divergence of the fields in the medium 
breaks the assumptions of the linear analysis and requires 
taking into account the dependence of the refractive index of 
medium on the light intensity. It was observed that simplest 
cubic Kerr-type nonlinearity does not lead to regularization 
of spectral singularities [117] but the saturable nonlinearity 
does [118]. 


5.3. Soliton scattering 

Collision of moving solitons with localized obstacles (de¬ 
fects) is the problem of fundamental significance in various 
nonlinear wave guiding systems. Recently in nonlinear op¬ 
tics a great deal of attention has been focused on the effect 
of soliton interaction with the complex PT-symmetric lo¬ 
calized potentials embedded into the self-focusing or self- 
defocusing medium. In these settings the PT-symmetric 
gain and loss regions are spatially localized, e.g, in the form 
of nonlinear PT-symmetric dimers, trimers, or, more gen¬ 
erally, oligomers embedded into nonlinear conservative lat¬ 
tice [44, 116,119,120]. Continuum versions of such models 
have been analyzed in the works [121-123]. In other cases 
soliton dynamics is investigated in the two-core systems 
such as an array of coupled PT-symmetric dimers [65,124] 
and its continuum version in the form of coupled planar 
waveguides with gain in one of them and balanced loss in 
another [59, 125]. Two-dimensional generalization of such 
models was studied in [126]. 


0.0018 



Figure 25 (a) Schematic plot of the waveguide array structure. 

The two PT-symmetric waveguides with balanced gain and loss 
are at the sites n = 0,1. The coupling constant between the 
modes of waveguides is C\ with the exception for the coupling 
constant C 2 between the active waveguides, (b) Total intensity 
of the localized mode excited by the incident soliton vs. p/p cr it 
for A = 0.2. (c,d) Intensity distribution in the waveguide array vs. 
the propagation distance in the process of soliton scattering for 
p = 1.5 = Pcrit/3 and the incident soliton amplitude (c) A = 0.2 
and (d) A = 0.5. Other parameters are C\ = 2, C 2 = 4, v = 0.5. 
Adapted from Refs. [119,120]. 


Soliton scattering by PT defects was first investigated 
in a nonlinear waveguide array [119,120]. A local inhomo¬ 
geneity was created by a pair of PT-symmetric waveguides 
with balanced gain and loss [see Fig. 25(a)]. Propagation of 
light in this system was modelled with the coupled discrete 
nonlinear Schrodinger equations 

i% + Cl(Vfc-l + Wn+ 1 ) + \Wn\ 2 Wn = 0, H ^ 0, 1, (50) 

m + /py/b + Ciy/_i+C 2 vo + |V'o| 2 V'b = 0, (51) 

m - ip V0+C 2 V'o + CiV'2 + |Vi|Vi = 0, (52) 


where the overdot stands for the derivative with respect 
to the propagation distance z, n is the waveguide number, 
i j/ n are the mode amplitudes at waveguides, p > 0(< 0) 
defines the rate of loss (gain) at zeroth and gain (loss) at first 
waveguides, and Cy2 are the coupling coefficients between 
the modes of waveguides. 

Far from the defect, i.e., for \n\ 1, Eq. (50) has an 

approximate soliton solution 

_ exp{z [v(ft — no) + (5 2 - v 2 )Ciz + a]} 
cosh [<5 (ft — no — 2Civz)] 


(53) 


which can be found by considering the continuum limit of 
Eq. (50) in the form of the integrable nonlinear Schrodinger 
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Figure 26 (a-b) Numerical results for the interaction of dark 

solitons with the PT-symmetric defect described by Eqs. (54,55). 
Shown is the distribution of | x P(x,z)|. (a) Reflection of slow dark 
soliton (v 0 = 0 . 1 ) and (b) the passage of the fast dark solitons 
(vo = 0.8). In all cases e = 0, xo = -10 y = 0.3. Adapted from 
Ref. [123]. (c,d) Scattering of gap solitons by PT defect according 
to Eqs. (57-59). Parameters are 77 = 5.86, ^/\rj \ = 0.02, E s = 
0.125, |v| = 0.05. Adapted from Ref. [122] 


equation. Here A, 8 =AyJ\ / (2C\), v, no, and a are param¬ 
eters defining the soliton amplitude, inverse width, velocity, 
initial position, and initial phase, respectively. 

It was found that solitons can be partially reflected and 
transmitted by the PT defect, and the total energy is gener¬ 
ally not conserved in the scattering process. In addition, the 
soliton can excite the mode localized at the PT-symmetric 
defect, and this process resonantly depends on the struc¬ 
ture parameters including the gain and loss, as illustrated in 
Fig. 25(b). Example of soliton scattering without excitation 
of localized mode is shown in Fig. 25(c), whereas under 
different parameters mode excitation is clearly visible, see 
Fig. 25(d). The soliton interaction with the large-amplitude 
localized mode can result in PT-symmetry breaking [120]. 

Interestingly, a discrete model describing a waveguide 
array with energy gain and loss that admits exact solitary 
wave solutions has been found [124]. Existence of the high- 
frequency and the low-frequency solitons was demonstrated. 
The numerical study of collisions between similar and dis¬ 
similar solitons has demonstrated that solitons of the same 
type interact elastically, while for the dissimilar ones the 
collision outcome significantly depends on the initial condi¬ 
tions. 

Karjanto et al. [123] have studied the scattering of bright 
and dark solitons on a strongly localized PT-symmetric 
defect in frame of the nonlinear Schrodinger equation for 


the local amplitude of the electromagnetic wave, T^z), 

PVz + + = [V(x) + iW(x )]'F, (54) 

with a complex potential: 

V (x) + iW (x) = e8(x) + iyS'(x), (55) 

where 8 and 8 ' denote the Dirac delta function and its deriva¬ 
tive, £ and y being real constants, positive or negative [123]. 
The nonlinear term in Eq. (54) represents defocusing non¬ 
linearity for g < 0. 

The collisions of incident dark soliton with the PT- 
symmetric defect were simulated, setting the initial con¬ 
ditions as 


¥(*, 0 ) = ^ cw 





vl(x - x 0 )\ + iv 0 


5 


(56) 

where xo is the initial soliton position and vo > 0 is its 
velocity. Shown in Figs. 26(a,b) are the examples of the 
interaction of the dark solitons with the defect for the cases 
of defocusing nonlinearity. In case (a) the soliton initial ve¬ 
locity is relatively small, and it is reflected by the defect. In 
case (b) the soliton velocity is relatively large, and it passes 
through the defect. An essential difference from previously 
studied interactions of solitons with defects in conservative 
systems is that, in spite of the gain-loss balance in the PT- 
symmetric dipole, the collisions change the norm of the 
solitons making the interaction dynamics more complex. 

Abdullaev et al. [122] studied scattering of gap solitons 
by PT defects in the framework of the following model, 


I'Pz+Vxx 

V 0 i(x) 

V d (x) 


<?| x P| 2x P + [V 0 i(x) + V’d(x)] x P, 


Vbcos(2x), 


rj+i^x / (x xq ) 2 \ 

2V )’ 


(57) 

(58) 

(59) 


that includes the periodic potential Voi(x) (optical lattice) 
with the amplitude Vo and the PT-symmetric localized defect 
Vd(v) that introduces gain and loss in the system. Here T] is 
the strength of the conservative part of the defect, coefficient 
£ stands for the gain-dissipation parameter, while the width 
of the defect is fixed to A = 5 in all numerical calculations. 
In the absence of the defect potential, Eq. (57) possesses 
families of exact gap soliton solutions. It is shown that by 
properly designing the amplitudes of real and imaginary 
parts of the PT defect it is possible to achieve a resonant 
transmission of the gap solitons through the defect. This 
phenomenon occurs for potential parameters that support 
localized modes inside the PT defect potential with the same 
energy and norm of the incoming soliton. When the imagi¬ 
nary amplitude of the PT defect is increased, the possibility 
of damping or amplification of transmitted or reflected gap 
soliton is demonstrated depending on the soliton incidence 
direction. Furthermore, unidirectional transmission or diode 
effect can be achieved as illustrated in Figs. 26(c,d). A pair 
of PT defects with oppositely positions gain and loss layers 
can then be used to trap a gap soliton. 
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Figure 27 (a) Schematic plot of the planar coupler having 

a defect of coupling exponentially localised around z = 0 and 
defined by Eq. (62). (b) Results for symmetric solitons (a = 1). 
Upper panels: field intensities with 77 = 0.15 interacting with defect 
characterised by K min = 1, kq = 2 and located at z = 0. Lower 
panel: respective evolution of the total energy flow So for / = 
1 (thick solid lines) and soliton amplitudes \q\\ and |^| (thin 
solid and dashed lines, respectively), (c) Same as in (b), but for 
antisymmetric solitons (<j = —1) and r\ = 0.5, / = 2.7, kq = 4. 
Adapted from Ref. [125]. 



Figure 28 (a) Modulation of gain and loss marked with shading 

along a waveguide, designed for irreversible coupling between two 
copropagating modes. Adopted from Ref. [139]. (b) PT-phase dia¬ 
gram showing the largest mode amplification vs. the normalised 
modulation frequency and gain/loss modulation strength. Modula¬ 
tion frequency ”2” corresponds to the regime in plot (a), where PT 
is in broken phase for arbitrary gain/loss. Adopted from Ref. [133]. 
(c,d) Irreversible coupling between modes of two waveguides with 
different cross-sections, and gain/loss modulation in the left wave¬ 
guide only. Light power evolution for input coupling to the (c) left 
and (d) right waveguide. Adopted from Ref. [140]. 


Soliton scattering on a defect in the two-core system 
schematically shown in Fig. 27(a) was studied by Bludov 
et al. [125]. The defect is introduced in the planar PT- 
symmetric coupled waveguides with gain (top waveguide) 
and loss (bottom waveguide) by local modification of the 
coupling constant in the vicinity of z = 0. Light propagation 
in this system can be described by the two coupled nonlinear 
Schrodinger equations 

r\ 

iq\, z = ~qi,xx + iyq 1 - K{z)q 2 - \qi \ q 1, (60) 

r\ 

iqi, z = -qi,xx - iyqi - k{ z)qi - ml qi, ( 61 ) 

K(z) = Ko-(Ko- K min )exp(— z 2 // 2 ). (62) 

The coupling coefficient k(z) attains the minimal value K m j n 
at z = 0 and tends to Kq at z — ±°°. The defect of cou¬ 
pling has the width /. Far from the defect (\z\ /) one 

has k(z) = Kb, and the model supports symmetric and an¬ 
tisymmetric types of soliton solutions as described above. 
Examples the of the soliton-defect interactions are presented 
for the cases of symmetric and antisymmetric solitons in 
Figs. 27(b) and (c), respectively. Upper panels give the held 
intensities, while lower panels present the respective evo¬ 
lution of the total energy flow So (thick solid lines) and 
soliton amplitudes \q\ | and |^| (thin solid and dashed lines, 
respectively). In Fig. 27(b), the soliton passes through a 
relatively short defect (/ = 1) transforming into a breather, 
which is characterized by the intensity oscillations between 


the two components: minimum (maximum) in one compo¬ 
nent corresponds to maximum (minimum) in the other. In 
Fig. 27(c), for l — 2.7 the splitting of the incident soliton 
in the two outward propagating pulses is observed. There 
also exist other scenarios of the soliton-defect interactions 
as discussed in Ref. [125]. 


6. Periodically Modulated Systems 

Periodic modulation is an effective tool for controlling tun¬ 
neling dynamics [127-131]. In recent years the PT sym¬ 
metry of the periodically modulated systems have attracted 
increasing attention. 


6.1. Modulated waveguide couplers 

Modulation of structure parameters along the propagation 
direction of coupled optical waveguides with gain and loss 
can open new opportunities for optical signal control in both 
linear [132-134] and nonlinear [50,62,135-138] regimes. 

Interestingly, the first studies of couplers with modulated 
loss and gain were carried out without a consideration of PT 
symmetry [139-142]. In these studies, it was shown that uni¬ 
directional energy transfer between co-propagating optical 
modes can be achieved in waveguides with longitudinally 
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modulated gain and loss regions. An example of a two-mode 
active coupler is shown in Fig. 28(a). Examples of unidi¬ 
rectional behaviour are presented in Figs. 28(c,d), where 
light gets localized and amplified in the right waveguide, 
irrespective of the input coupling port. Such irreversible cou¬ 
pling between the co-propagating modes cannot be achieved 
in conservative structures. This phenomenon relies on the 
breaking of the time-reversal symmetry through the modula¬ 
tion of the complex optical refractive index, being analogous 
to unidirectional scattering by PT-symmetric Bragg gratings 
discussed in Sec. 5.1. 

Light evolution in modulated two-mode couplers can be 
described within a framework of coupled-mode theory, 

d_ (a\\ = ( iyi(z)-Pi(z) -C(z) \ / a A (63) 

dz\a 2 ) V ~C(z) in(z)-P2(z)J\a 2 J 

where z is the propagation distance along the waveguide, 
cij are the mode amplitudes, C is the coupling coefficient, 
Yj are the gain/loss coefficients, are the real propagation 
constants, and j = 1,2 is the mode number. Such type of 
models were found to possess PT symmetry, under certain 
restrictions on the modulation profiles [132-134]. 

The formulation of Eq. (63) is not unique, and it de¬ 
pends on the choice of the basis modes. For the structure 
presented in Fig. 28(a), if the mode amplitudes are chosen 
as a\ = < 2 + + a _ and <22 = a + — a_, where a± correspond 
to symmetric and antisymmetric supermodes in the ab¬ 
sence of gain or loss, we have Yi(z) = — Y\ (z), C(z) = C(0), 
Pj(z) = /3o- Such model is PT-symmetric in case of even 
modulation, Yj(z) = Yj{~z). PT diagram for harmonic mod¬ 
ulation, yi = — 72 ~ cos(o)z) was calculated in Ref. [133], 
see Fig. 28(b). In PT-symmetric regime, average amplifi¬ 
cation is zero, however in PT-broken case, one mode ex¬ 
periences strong amplification. Boundaries between these 
regimes depend nontrivially both on the modulation am¬ 
plitude and frequency. In case of resonant modulation am¬ 
plitude matching the beating period between the modes, 
which corresponds to the normalized value “ 2 ” along the 
vertical axis in Fig. 28(b), PT mode symmetry is broken for 
arbitrarily small values of gain/loss modulation. It is in this 
regime that unidirectional mode coupling is always realized, 
as illustrated in Figs. 28(a,c,d). This is a key difference from 
unmodulated linear PT couplers, where PT breaking occurs 
only above a certain threshold of gain/loss. 

Nonlinear optical interactions can have a nontrivial ef¬ 
fect on the mode evolution. It was shown that local Kerr- 
type nonlinearity can preserve or change the stability prop¬ 
erties. Of particular interest is the regime when the sys¬ 
tem is unstable in the linear regime, but becomes stable 
due to conservative nonlinear self-action, as illustrated in 
Fig. 29(a). Such stabilising effect can enable design of laser 
sources, similar to the case of actively coupled PT wave¬ 
guides with fixed parameters discussed in Sec. 2.3, however 
periodic modulation can be realized in a simple way whily 
allowing fine control over the emerging nonlinear dynam¬ 
ics, which can be periodic or aperiodic. Indeed, chaotic 
evolution was predicted for a coupler (dimer) with period¬ 
ically modulated gain and loss [138]. This is illustrated in 


Figure 29 Dynamics in nonlinear modulated PT structures. 

(a) Evolution of intensities in a PT dimer with periodically modu¬ 
lated coupling. Parameters correspond to nonlinear stabilisation 
in a linearly unstable system, adopted from Fig. 8 in Ref. [137]. 

(b) Stroboscopic Poincare map for a dimer with periodically mod¬ 
ulated gain and loss. Adopted from Ref. [138]. (c,d) Dynamics 
in a trimer demonstrating (c) stable and (d) unstable evolution. 
Adopted from Ref. [136]. 


Fig. 29(b), where the phase plane is presented for variables 
(V',V')> with v)(z) = | ^ 2 (z) P — \&i (z) | 2 . Periodic modula¬ 
tion can be further applied to multiple coupled nonlinear 
waveguides, or oligomers [50,136]. Examples of stable pe¬ 
riodic and unstable dynamics in a nonlinear trimer with 
periodically modulated gain and loss are show in Figs. 29(c) 
and (d), respectively. Furthermore, periodic modulation can 
stabilize temporal solitons in coupled waveguides, in the 
supersymmetric regime when loss and gain are at the linear 
PT breaking threshold [135]. 

When system parameters are modulated such that un¬ 
derlying system is not PT-symmetric at any modulation 
cross-section, the overall dynamics can still demonstrate 
the features of pseudo-PT symmetry [143]. The concept of 
pseudo-PT symmetry corresponds to the PT-symmetry in 
the effective system, which manifests as quasi-stationary 
propagation. Such regime can be achieved under high- 
frequency periodic modulation, when the modulated system 
can be mapped into an effectively unmodulated one with 
rescaled parameters. Following Ref. [143], consider a cou¬ 
pler described by Eq. (63) with balanced gain and loss, 
Yi(z) = 7/2 and 72 (z) = - 7 /2, fixed coupling C(z) = — v, 
and modulated propagation constants Pi(z) = —S(z)/2, 
(z) = S(z)/ 2. A biharmonic modulation was analyzed, 
S(z) = — -A[sin(fi)z) + /sin(2o)z + </>)]. For </) ^ 0 , 7 r, the 
corresponding Hamiltonian is not PT-symmetric. Never¬ 
theless, the averaged system obtained after averaging the 
high-frequency modulations features PT symmetry, and the 
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Figure 30 Quasi-PT symmetry in linear structures, (a-c) Adopted 
from Ref. [143]. (a) Schematic of a coupler with symmetric gain 
and loss, but with asymmetric modulation of the propagation con¬ 
stants. (b) The mode growth rate vs. the normalized gain/loss am¬ 
plitude (horizontal axis) and the modulation strength (vertical axis), 
for fast modulations. White line marks the boundary of quasi-PT 
transition, (c) Example of almost periodic evolution with very small 
growth rate in quasi-PT symmetric regime, (d) Phase diagram 
analogous to Hofstadters butterfly spectrum for a PT-symmetric 
Hamiltonian with aperiodic driving, adopted from Ref. [134]. 


quasi-energies for such effective system are found as 
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Figure 31 (a) Schematic of a waveguide array with modulated 

gain and loss according to Eq. (66). (b) Phase diagram of PT mode 
symmetry. The inset shows the real (upper panel) and imaginary 
part (lower panel) of the mode quasi-energies, (c) Beam intensity 
profiles for (c) conventional propagation across the lattice in a 
conservative array and (b) hyper-ballistic regime in a modulated 
lattice close to PT threshold. Adopted from Ref. [144]. 


agation and spreading through the lattice [144]. Detailed 
theoretical analysis was performed for a lattice with periodic 
modulation of gain and loss, which is out-of-phase between 
the neighbouring waveguides as illustrated in Fig. 31(a). 
Linear beam evolution in such arrays is described by the 
coupled-mode equations, 


£ = ±\j\yf\ -[y/(2/)] 2 , 

with the rescaled coupling strength 


(64) 


da 

i~r~ = -K(a n ~ i + a n+l ) + i(-l) n A(t)a n , (66) 

dz 


(s)-'"(to) exp( '""'' ) - <65) 

The two quasi-energies are real if y < 2|/|, and they are 
complex if y > 2\J\ with the critical point y c = 2\J\. There 
exists a spontaneous PT-symmetry-breaking transition in 
the effective system when the imaginary part Im(e) changes 
from zero to nonzero, as illustrated in Fig. 30(b). It appears 
surprising that the quasi-energies can be real even when 
the modulated system (63) is non-PT-symmetric, however 
the numerical analysis reveals that the eigenvalues of the 
original modulated structure have a very small but nonzero 
imaginary part. Accordingly, at intermediate propagation 
distances the dynamics can be almost exactly periodic in 
the quasi-PT symmetric regime, as show in Fig. 30(c). 

Further interesting possibilities can arise in case of ape¬ 
riodic driving. It was found that a phase diagram analogous 
to Hofstadters butterfly spectrum can emerge in PT Hamil¬ 
tonians [134], see Fig. 30(d). The horizontal axis in the plot 
corresponds to the ration of two driving frequencies, and 
the vertical axis characterizes detuning between the modes. 

6.2. Modulated lattices 

Arrays of coupled modulated waveguides with gain and loss 
open new opportunities to control the rate of beam prop- 


where n is the waveguide number, a n are the mode am¬ 
plitudes, t is the propagation coordinate, K is the cou¬ 
pling between the neighbouring waveguides, and A (t) = 
Ao square (cot + 0) is a square modulation profile of gain 
and loss. Solutions of this system can be found analytically 
in Fourier space, characterized by the transverse wavenum¬ 
ber q. For a range of modulation parameters indicated with 
shading in Fig. 31(b), the eigenmodes with all q have unbro¬ 
ken PT symmetry. Outside of this parameter region, some 
modes have broken PT symmetry. Characteristic dispersion 
dependencies corresponding to these two regimes are shown 
as insets in Fig. 31(b). The imaginary part at the lower insets 
represents the mode gain or loss, and we see that it is zero for 
all modes in the unbroken PT regime. The real part shown in 
the upper insets defines the propagation angle (proportional 
to the first derivative) and strength of diffraction (propor¬ 
tional to the second derivative). It is seen that these propaga¬ 
tion characteristics can be very strongly affected due to PT 
transition. A comparison of beam propagation in a conser¬ 
vative array and array with gain/loss modulation parameters 
chosen close to PT threshold are presented in Figs. 31(c) 
and (d), respectively. The input beam is inclined such that 
the phase gradient between the neighbouring waveguides is 
Ti/ 4, which corresponds to the fastest propagation across a 
conventional conservative lattice, limited by the coupling 
rate K [Fig. 31(c)]. However in a PT lattice with modulated 
gain and loss, hyper-ballistic transport can be realised with 


J=v £ (i i)~ m J-2m 

m=—oo 
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Figure 32 (a) Schematic of a PT waveguide array with bending 

loss in every second waveguide, (b) The eigenvalue spectrum 
vs. the transverse lattice wavenumber q. (c) Experimental fluo¬ 
rescence microscopy image of the light beam as propagating 
through the lattice, (d) The variance of the beam vs. the propaga¬ 
tion distance z (red line, experimental data; black line, simulation). 
Adopted from Ref. [145]. 


much higher propagation velocity [Fig. 31(c)]. Additionally, 
the beam spreading due to diffraction is strongly reduced. 

Experimental realisation of PT lattices based on fs laser- 
written modulated conservative waveguides was reported 
in Ref. [145]. As shown in Fig. 32, every second wave¬ 
guide was periodically curved, which resulted in bending 
losses. In the regime of high-frequency modulation, the 
loss can be considered as effectively constant along the 
waveguides. The corresponding spectrum of eigenmodes 
is presented in Fig. 32(b). Eigenmodes in different regions 
of the spectrum will experience different losses, depend¬ 
ing on the value of the imaginary part of their propagation 
constant. As a consequence, in both bands modes in the 
centre of the spectrum (where q ~ 0) experience an interme¬ 
diate loss, which is close to the average loss in the system. 
At the edge of the spectrum (around q ~ 7r), the modes in 
the upper band suffer from a loss that is much higher than 
the systemss average loss, whereas in the lower band the 
modes experience much less loss. Hence, because of decay, 
the modes in the upper band at the edge of the spectrum 
will disappear after a relatively short propagation distance 
z, whereas the modes in the centre of the spectrum will 
disappear somewhat later. Therefore, at long propagation 
distances, the spectrum will be getting narrower, and only a 
part of the spectrum will contribute to transport. Experimen¬ 
tal image of the beam propagation is presented in Fig. 32(c), 
and Fig. 32(d) presents the beam variance characterising 
its width, o 2 (z) = Hn n 2 \ a n \ 2 / Hn \ a n \ 2 ■ It is observed that 
initially, the wave exhibits ballistic spreading with o 2 ~ z 2 , 
as all eigenmodes still contribute to the transverse, transport. 
However after z C rit — 2cm, transition from ballistic to dif¬ 
fusive transport is visible, with slower beam spreading as 
a 2 ~ z. Such coexistence of ballistic and diffusive transport 
in a static, ordered system, could only happen due to non- 
Hermitically. More recently, topological phase transition 
was reported in such waveguide arrays [146]. 
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Figure 33 (a) Schematic diagram of the temporal fibre networks. 
Two coupled fibre loops periodically switching between gain and 
loss as used in the experiment. Pulses are delayed or advanced 
owing to a length difference AL between the short and long loops, 
(b) The evolution of the pulses in the temporal networks. Passages 
through the short and long loops are indicated, (c) The equivalent 
PT-symmetric spatial mesh lattices. Gain (red) and loss (blue) 
channels are arranged anti-symmetrically and are periodically 
coupled. The real part of the potential is introduced by using phase 
modulation with ±<po. From Ref. [148]. (d) and (f): Simulation and 
experiment results in the nonlinear regime above PT threshold, (e) 
and (g): Simulation and experiment results in the nonlinear regime 
below PT threshold. Optical soliton and internal oscillations are 
observed in the experiment. From Ref. [149]. 


The unconventional properties of modulated PT lattices 
based on waveguide arrays suggest that such systems would 
demonstrate novel features in the nonlinear regime, and 
we anticipate theoretical and experimental studies targeting 
these topics in the near future. In particulatr, due to the 
flexibility of waveguide array writing with fs lasers [147], 
this platform could enable the observation of linear and 
nonlinear effects in modulated PT structures. 


6.3. Mesh lattices 

Regensburger et al. [148] demonstrated the first experimen¬ 
tal realization of a large-scale PT-symmetric lattice, which 
was a new kind of optical PT synthetic devices. They re¬ 
ported the experimental observation of light transport in 
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large-scale temporal lattices which are parity-time symmet¬ 
ric. In addition, the periodic structures respecting PT symme¬ 
try can act as unidirectional invisible media when operated 
near their exceptional points. 

To realize the specially artificial structures and study 
the light transportation, they have used an elegant experi¬ 
mental arrangement that operates in the temporal domain, 
see Fig. 33 (a) and (b). In most cases, the refractive index 
is a complex quantity. The real part of the refractive index 
affects the velocity of the propagating light, while the imagi¬ 
nary part can lead to the amplification or absorbtion of light 
within a material. By temporally modulating the real and 
imaginary part of the refractive index, they were able to 
realize the PT-symmetric optical network in the temporal 
domain, whose system is analogous to a spatially periodic 
mesh network with gain and loss. Each node of the spatial 
network corresponds to a specific ’time slot’ of the tem¬ 
poral lattice, see Fig. 33 (c). A sequence of light pulses 
are injected into two connected optical fiber loops that are 
designed to exhibit PT symmetry. The anti-symmetric imagi¬ 
nary part of refractive index profile is attained by alternating 
gain and loss in the two loops via using optical amplifiers 
and amplitude modulators. While the even, real component 
of the index profile is introduced by using phase modulators. 

The observation of the unusual unidirectional invisibil¬ 
ity in their system is also interesting. Since the left-right 
symmetry is broken in this structure and propagation is no 
longer invariant when gain and loss are exchanged in time. 
More specifically, light propagating in such a system can 
experience reduced or enhanced reflections dependent on 
the direction of propagation while it has little effect on the 
transmitted light. The system can become totally invisible 
when light traverses it from one side, whereas it can still be 
seen when it is illuminated from the other. 

In Ref. [150], they investigated a class of optical mesh 
periodic structures that are discretized in both the transverse 
and longitudinal directions, which is equivalent to the above 
ones in the temporal domain. The mesh arrangements are 
composed of an array of waveguides with each one being 
discretely and periodically coupled to its adjacent neigh¬ 
bors. The practical appealing of this type of lattice is the 
physical separation between the coupling and amplification 
within each building block. And the band structure of these 
systems have been systematically analyzed and their dis¬ 
persion relation is analytically obtained. And the optical 
dynamics including the unidirectional invisibility in these 
PT-symmetric mesh lattices are also examined. 

By introducing nonlinearity, the observation of stable 
optical discrete solitons in the PT-symmetric mesh lattices 
system has been experimentally demonstrated [149]. The 
structure of lattices is the same as shown in Fig. 33 (c), 
which is globally PT-symmetric. In the linear regime, when 
a single position is excited, the wavepacket spreading ex¬ 
periences a great amplification in the broken PT-symmetry 
regime and remains neutral below the PT threshold. Inter¬ 
estingly, when the power of single initial pulse is raised, the 
nonlinearity comes into play and the discrete solitons may 
appear in the system. Above the PT threshold, the nonlinear¬ 
ity cannot lead to stable solitary waves and the wavepacket 


still spreads exponentially, see Fig. 33 (d) and (f). However, 
below the PT threshold, the nonlinearity can compensate 
the spreading and the optical soliton appear, see Fig. 33 (e) 
and (g). It has been demonstrated that these solitons belong 
to a continuous family of solitons which is not typically the 
case for dissipative solitary waves and this class of discrete 
solitons can avoid instabilities. Based on such mesh lattices 
system, the behavior of nonlinear waves provides a possibil¬ 
ity for realizing saturable absorbers which are widely used 
in Q-switched laser cavities and ultra-short optical pulse 
arrangements. 

These results, obtained for both temporal and spatial 
modulated systems, represent an important step towards 
applications of the parity-time symmetry concept to a new 
generation of multi-functional optical devices and networks. 

7. Conclusion and outlook 

Future technologies will demand a substantial increase in a 
density of photonic integration and energy efficiency far sur¬ 
passing that of bulk optical components and modern silicon 
photonics. Such advances can be achieved only by embed¬ 
ding the photonic functionalities at materials level taking 
advantages of sub wavelength nanophotonics and metade¬ 
vices. Many metadevices based on plasmonic technologies 
would require an engineered loss compensation that, as we 
expect, can be achieved in a clever way through the con¬ 
cepts of non-Hermitian dynamics and the PT symmetric 
photonic systems. It is believed that the PT symmetric pho¬ 
tonic systems can offer promising solutions of a range of 
important problems providing an efficient and smart loss 
compensation. Photonic PT systems feature alternating and 
symmetrically distributed regions of gain and loss, and in 
such systems losses are as important as gain so that, depend¬ 
ing on the structural parameters, gain can either compensate 
losses or amplify optical pulses. In addition, tunable pho¬ 
tonic PT-symmetric systems can add many new exciting 
functionalities being potentially able to find even broader 
applications beyond the field of photonics, e.g. in the field 
of Bose-Einstein condensates [151]. Also, the soliton scat¬ 
tering on PT-symmetric defects has been studied for otehr 
types of models, for example for the Klein-Gordon sys¬ 
tem [152,153]. 

As an example of the most recent developments we 
mention a very general study of constant-intensity waves in 
non-Hermitian potentials. It is well known that any finite- 
amplitude wave propagating in a nonlinear system under¬ 
goes modulational instability in any Hermitian potential 
with self-focusing nonlinearity. Markis et al. [154] demon¬ 
strated that this fundamental restriction is conveniently lifted 
when working with non-Hermitian potentials. In particular, 
they presented a whole class of waves that have constant 
intensity in the presence of linear and nonlinear inhomoge¬ 
neous media with gain and loss. These results suggest new 
directions for the experiments on nonlinear non-Hermitian 
scattering. 

Many recent developments on non-Hermitian and PT 
symmetric systems occur beyond the field of photonics, and 
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one of the examples is exciton-polariton condensate which 
is composed of hybrid light-matter quasiparticles formed 
by strongly interacting photons and excitons (electron- 
hole pairs) in semiconductor microcavities. The exciton- 
polaritons always exist in a balanced potential landscape 
of gain and loss and, as was shown recently [155], the non- 
Hermitian nature of the system can modify dramatically 
the structure of modes and spectral degeneracies in exciton- 
polariton systems, and, therefore, will affect their quantum 
transport, localisation, and dynamical properties. These find¬ 
ings pave the way for studies of non-Hermitian quantum 
dynamics of exciton-polaritons, which can uncover novel 
operating principles for polariton-based devices. 

As discussed above, PT symmetric systems demon¬ 
strate many nontrivial non-conservative wave interaction 
and phase transitions, which can be employed for signal 
filtering and switching, opening new prospects for active 
control of light. In this review, we have discussed only some 
problems involving nonlinear PT-symmetric photonic sys¬ 
tems with an intensity-depend refractive index. Nonlinearity 
in such PT symmetric systems provides a basis for many 
effects such as formation of localized modes, nonlinearly 
induced PT-symmetry breaking, and all-optical switching of 
signals. Nonlinear PT-symmetric systems can serve as pow¬ 
erful building blocks for the development of novel photonic 
devices targeting an active light control. 

As follows from the results summarized above, the study 
of photonic systems with the PT symmetry is an active re¬ 
search area which is expanding very rapidly leaving the 
boundaries outlined by the photonic systems and penetrat¬ 
ing to other fields such as the physics of Bose-Einstein con¬ 
densates. However, the photonics is a cornerstone of those 
research due to rapidly expanding experimental demonstra¬ 
tions. Here, we have reviewed only a few of many intrigu¬ 
ing properties and applications of PT-symmetric photonics 
systems being concentrated on the less studied nonlinear 
properties, but the field is quickly growing and many new 
ideas and demonstrations appear almost every week. 
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